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We compute the lattice 1-site probabilities, on the upper half-plane, of the four height 
variables in the two-dimensional Abclian sandpile model. We find their exact scaling form 
when the insertion point is far from the boundary, and when the boundary is either open 
or closed. Comparing with the predictions of a logarithmic conformal theory with central 
charge c — —2, we find a full compatibility with the following field assignments: the heights 
2, 3 and 4 behave like (an unusual realization of) the logarithmic partner of a primary 
field with scaling dimension 2, the primary field itself being associated with the height 1 
variable. Finite size corrections are also computed and successfully compared with numerical 
simulations. Relying on these field assignments, we formulate a conjecture for the scaling 
form of the lattice 2-point correlations of the height variables on the plane, which remain as 
yet unknown. The way conformal invariance is realized in this system points to a local field 
theory with c = — 2 which is different from the triplet theory. 

PACS numbers: 05.65. -Hb,45.70.-n 

I. INTRODUCTION 

Sandpile models are open dynamical systems which generically show a very rich spectrum of 
critical behaviours, both spatially and temporally. Along with many other models, they are complex 
critical models. Because of their non-linear dynamics and their non-local features, a fixed external 
driving can trigger events whose scales follow critical distributions. The references [1-3] review the 
ideas and models involved. 

One of the simplest and yet most challenging models is the two-dimensional model originally 
proposed by Bak, Tang and Wiesenfeld [4], now referred to as the BTW model, or the unoriented 
Abelian sandpile model (ASM) after it was shown by Dhar to possess an Abelian group [5]. It is 
this model we consider in this article. 

The most natural formulation of the ASM is in terms of discrete height variables, taking the 
four integral values 1 to 4 and located at the sites of a grid. The configuration of the sandpile, 
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formed by the values of the heights, evolves under a stochastic dynamics, which eventually reaches 
a stationary regime, where the configurations are weighted according to an invariant measure. 
Numerical simulations showed that this model has critical properties [4]. 

The first analytic treatment of the ASM was carried out by Dhar [5], where general and im- 
portant features of the ASM dynamics were obtained. That work paved the way for numerous 
exact results. Among these, one can mention the results on correlations of height 1 variables and 
related clusters (the so-called weakly allowed sub-configurations) [6-10], on 1-site probabilities of 
height variables on the whole plane [11], on boundary correlations of height variables [12-14], on 
avalanche distributions [15-19], on specific effects of boundary conditions [20-22], on finite size 
corrections [20, 23], and on the insertion of dissipation [8, 13, 21, 22]. Moreover, a number of these 
references analyze their results in the light of a possible conformal field description. All available 
results are compatible with such a description, and point to a specific logarithmic conformal field 
theory, having central charge c = — 2. 

Quite recently, the 1-site probability for the height 2 variable on the upper half-plane, in fact a 
disguised 2-point correlator, has been computed exactly in the scaling limit, the results having been 
reported in [24]. It was proved, as had been suspected before, that the height 2 variable behaves 
logarithmically, unlike the height 1: the height 2 correlations contain logarithms, while those for 
the height 1 variables are rational functions of the distances. Moreover, numerical simulations 
suggested that the height 3 and 4 variables have the same behaviour and scale in the same way 
as the height 2. Finally, the logarithmic conformal field to which these variables converge in the 
scaling limit was identified, and then used to make predictions for the 2-site bulk correlations of 
height variables, which, to this date, have not been computed on the lattice (except for height 1). 
In the logarithmic conformal field theory terminology, the scaling limit of the height 2 (and 3 and 
4) is a logarithmic field, which is the partner of the field to which the height 1 scales. 

Our purpose in this article is two-fold. First, we want to give the details of the analysis reported 
in [24]. It involves a lattice part, where the probability of interest is computed in the ASM, and 
a conformal theory part, where the relevant 4-point correlation function with the intended scaling 
field is calculated and then compared with the lattice results (and simulations). The details on 
both, too long and therefore omitted in [24], will be given here. Second, we have extended the 
analysis to include the height 3 and 4 variables. Our new exact results confirm the behaviour 
conjectured in [24]. Since the height 3 and 4 variables each scale like linear combinations of the 
height 1 and 2 variables, any multisite height probability can be reduced, in the scaling limit, to 
multisite probabilities involving heights 1 and 2 only. These results should close the long chapter 
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devoted to the height probabihties in the ASM and their description by a field theory . 

That the critical properties of the ASM can indeed be described by a conformal field theory re- 
mains for us a major motivation, for two reasons. First, the description by an underlying conformal 
field theory can be considered as an actual and explicit solution of the lattice model. Second, the 
ASM would provide the first instance where the correspondence between a logarithmic conformal 
field theory and a lattice realization is most clearly exposed. We believe that such a correspon- 
dence, so common for usual equilibrium lattice models, should inevitably clarify tricky issues on 
both sides. 

The article is organized as follows. In Section II, we define the sandpile model we consider, 
review its basic features and recall the elements which prove useful for what follows. 

In Sections III and IV, we review the graph theoretic tools used for the derivation of the 1-site 
height probabilities on the plane, in a way that allows the formalism to be applied to the upper 
half-plane (UHP). The calculation presented in Section III for the height 2 is identical to that of 
Priezzhev [11], but our derivation for the height 3 in Section IV is different. The method we use 
leads to the same results as Priezzhev's but has an interesting corollary, since it leads to an exact 
relation between the probabilities P2 and P3, a relation unnoticed so far. Moreover, based on very 
strong numerical evidence, we conjecture exact and simple formulas for P2 and P3, similar to that 
for Pi , namely, simple polynomials in ^ . Although we could not provide a proof for them, we have 
checked that these formulas are accurate to one part in 10^^. 

Section V applies the formalism developed in the previous two sections to the upper half-plane, 
with two homogeneous boundary conditions on the real axis: open and closed. The resulting 
formulas are then collected in Section VI: the four probabilities P^^m), Pf(rn), P2^(m) and 
P3^(m), where m is the distance of the reference site to the boundary, are explicit but complicated 
expressions, involving multiple integrals and summations. 

In Section VII, we proceed to the asymptotic analysis of the integral expressions just derived. 
These expressions have expansions in inverse powers of m, with coefficients which are polynomials 
in logm. Anticipating that the scaling fields corresponding to the height variables have dimension 2, 
we restrict the analysis to the terms of order 2 in m~^. The dominant terms in the probabilities are 
of order m~^, as expected, and are computed exactly. The final results are summarized in Section 
VIII. The reader interested in the results, but not in the technical details of the calculations, can 



^ With the usual reservations: we have identified the fields describing the height variables provided that indeed such 
an identification exists at all, i.e. is compatible with every single lattice joint probability. 
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safely go straight to Section VIII. 

Section IX starts the conformal field theoretic side of our work. In order to assess the effect of a 
change of boundary condition on the 1-site probabilities, we first compute the correlation function 
{n{zi)fi{z2)ipiz, z)) on the UHP, where the /i fields change the boundary condition at two points 
zi and Z2 on the real axis, and where the non-chiral field il^{z, z) is supposed to describe the scaling 
limit of the height variables (2, 3, or 4). The above 3-point field correlator enables us to relate to 
each other the two UHP 1-point functions (■0) in front of a boundary which is either fully open or 
fully closed. When if) is taken to be the logarithmic partner of a dimension 2 primary field ^, this 
relation exactly reproduces the way the probabilities P°^(m) and Pf{m) are related on the lattice. 

The identification with -0 is further confirmed in Section X, where the finite size corrections 
predicted by the conformal theory are compared to the results of numerical simulations of the 
lattice model, showing an excellent agreement. 

Section XI reconsiders the field correlation functions established before from a purely conformal 
field theory point of view. The local logarithmic conformal theory with c = —2, believed to be 
relevant here, has a well-known realization in terms of free symplectic fermions 6, [25-27]. The 
primary field (j) describing the height 1 variable in the ASM has been well studied and is known 
to be expressible as a local field in the 9, 9 theory. We show that the field ip, corresponding to the 
heights 2, 3 and 4, has no local realization in the 6, 9 theory. It follows that the well-studied triplet 
theory with c = — 2 [26] is not the appropriate logarithmic conformal field theory to describe the 
sandpile model. 

Finally, relying on the field identifications obtained before, we discuss in Section XII the height 
joint probabilities on the plane, and formulate conjectures for the 2-site lattice probabilities on the 
plane, unknown so far. These conjectures compare well with the results obtained from simulations. 

II. THE ABELIAN SANDPILE MODEL 

We consider a finite, two-dimensional grid C and attach to each site i a random variable hi, 
representing the height of sand at that site and taking integer values 1, 2, 3, ... A configuration C 
of the sandpile is the set of values {hi} for all sites. The probability distribution over the height 
variables or equivalently, over the configurations, evolves dynamically from some initial distribution 
Po(C). 

The discrete dynamics, itself stochastic, is completely defined in terms of a toppling matrix A. 
This matrix has integral entries, strictly positive on the diagonal, negative or zero off the diagonal, 
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and moreover satisfies the condition that its row sums Ylj ^ij non-negative. A configuration 
is called stable if its height values satisfy hi < An for all i. 

At time t, the dynamics takes the stable configuration Ct into a new stable configuration Ct+i 
as follows. To the configuration Ct, we first add a grain of sand at a random site i by setting 
hi ^ hi + 1 (we assume that the site i is chosen randomly with a uniform distribution on C). This 
new configuration, if stable, defines Ct+i- If it is not stable, the site i which has been seeded has 
now a height hi = An + 1 > An, and, as consequence, topples: it loses An grains of sand, the 
other sites j each receive — Ay grains, and Ay grains are dissipated. In other words, when the 
site i topples, the heights change according to 

hj ^ hj - Aij, yj e L. (2.1) 

After the site i has toppled, other sites can be unstable (have heights hj > Ajj), in which 
case they too topple according to the same toppling rule (2.1). Once all unstable sites have been 
toppled, one is left with a stable configuration, Ct+i- One can show [5] that provided there are 
dissipative sites in C (namely, sites k for which J2j ^kj > 0), this dynamics is well-defined: the 
order in which the unstable sites are toppled does not matter, and the new stable configuration 
Ct+i is reached after a finite number of topplings. 

In the usual BTW model, which we consider here, the grid £ is a rectangular portion of the 
discrete plane Z^, and A is the discrete Laplacian on C with appropriate boundary conditions on 
dC Each boundary site k can be either open or closed, depending on whether it is dissipative or 
not. Explicitly, we take the toppling matrix to be 

Zi if i = j is either a bulk site or a closed boundary site, 
4 a i = j is an open boundary site, 

(2.2) 

—1 if i,j are nearest neighbours. 



Aij = < 



otherwise, 

where Zi is the coordination number of i in C, equal to 4 for a bulk site, to 3 for a boundary site 
and to 2 for a corner site. The toppling rule (2.1) then says that when a bulk or a closed boundary 
site topples, the sand grains it loses are all transferred to the nearest neighbours, and when an 
open boundary site topples, it loses 4 grains of sand, Zi of them going to the nearest neighbours, 
and the others falling off the sandpile. In the latter case, it is convenient to think of the sand grains 
falling off the pile as being transferred to a sink site, connected to all dissipative sites. Bulk and 
closed boundary sites are conservative {J2j A^j = 0), while open boundary sites are dissipative 
(S^.A,,>0). 
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Because it involves adding a sand grain at a random site, the dynamics is stochastic and can be 
written in terms of a probabihty transition matrix p. The matrix p gives rise to a master equation, 
Pf+i(C) = J2c' Pc,C' PtiC), which specifies the time evolution of the probability measure on the 
configurations. 

The long-time behaviour of the sandpile is described by the time-invariant probability measures. 
For the BTW model as defined above, Dhar has shown that there is a unique time-invariant measure 

[5]. This measure assigns all stable configurations their probability of occurrence, when the 
dynamics has been applied for a sufficiently long time so that the system has set in the stationary 
regime. The thermodynamic limit of is what we want to compare with a conformal field 
theoretic measure. 

The number of stable configurations is Yli but only a small fraction of them keep reappearing 
under the dynamical evolution. The transient configurations are not in the repeated image of the 
dynamics, and occur only a finite number of times. As a consequence, they all have a zero measure 
with respect to P^. A simple example of a transient configuration is /ij = 1 for all i. The non- 
transient configurations are called recurrent and asymptotically occur with a non-zero probability. 

For the BTW model, the recurrent configurations all occur with equal probability under the 
dynamics [5]. The invariant measure P^ is thus a simple uniform distribution, but its support, 
the set of all recurrent configurations, is not simple, which explains why exact calculations in the 
ASM are notoriously difficult. On general grounds, it can be shown that the number of recurrent 
configurations is equal to = det A [5]. For the particular model discussed here, this number is 
roughly equal to 3.2ll^l (up to surface terms), an exponentially small fraction of the 4l^l stable 
configurations. 

Transient and recurrent configurations can be distinguished by a criterion based on forbidden 
sub-configurations (FSC's): a subset of sites /C C £ is said to form a forbidden sub-configuration 
if hi < — J2j(^ic\{i} '^ji ~ #{neighbours of i in /C}, for all sites i of /C. For instance, two neighbour 
sites with height 1 form an FSC. Then a configuration is recurrent if and only if it contains no 

FSC [23]. 

A practical way to test whether a configuration is recurrent or transient is to use the burning al- 
gorithm [5] . At time 0, the sink site is the only site to be burnt: we define /Co = >C to be the set of un- 
burnt sites at time 0. The sites i of JCq such that hi > — J2jeKo\{i} ~ #{neighbours of i in ACo} 
are burnable at time 1 (they can only be dissipative sites, namely, open boundary sites). So we 
burn them, obtaining a set /Ci C /Cq of unburnt sites at time 1. Then the sites of /Ci which are 
burnable at time 2, i.e. those satisfying hi > —J2j^!Ci\{i} ~ #{neighbours of i in /Ci}, are 
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burnt. This leaves a set /C2 C /Ci of unburnt sites at time 2. This burning process is carried 
on until no more sites are burnable, which means that /Ct+i = K,t for a certain T. Then the 
configuration is recurrent if and only if all sites of C have been burnt (/Ct = 0)- Otherwise Kt is 
an FSC, and the configuration is transient. 

The burning algorithm allows the definition of a unique rooted spanning tree on a graph £*, 
from the path followed by the fire in the lattice [23] . The graph L* has the sites of L and the sink 
as vertices, and has links defined by A: an off-diagonal entry Ajj = — n means there are n bonds 
connecting the sites i and j, and each site i is connected to the sink by a number of bonds equal 
to Y^j^c — ^' namely the number of grains dissipated when i topples. At time 0, the sink is 
the only burnt site and forms the root of the tree. In the next steps, the fire propagates from the 
sink to those sites which are burnable at time 1, then from the sites which have been burnt at 
time 1 to those which are burnable at time 2, and so on. If a site burns at time t, it catches fire 
from one among its neighbours that were burnt at time t — 1 (or from the sink site at time 1). In 
case there are more than one of these, a fixed ordering prescription is used to decide along which 
bond the fire actually propagates [23]. The collection of all bonds forming the fire path defines a 
spanning tree, rooted in the sink, and growing towards the interior of the lattice JC. Conversely, to 
every spanning tree — from Kirchhoff 's theorem, their number also equals Af = det A — one can 
associate a unique recurrent configuration if the ordering prescription mentioned above is known. 

Since there is a one-to-one mapping between the recurrent configurations and the spanning 
trees, the latter provide alternative variables for the description of the ASM. So we trade the 
random height variables for random spanning trees, on which the invariant measure is uniform. 
Although the local height variables look more natural, they are strongly correlated over the whole 
lattice, since their values have to make up a recurrent configuration. Their long-range correlations 
are encoded in the global structure of the spanning trees. Thus height variables are local but 
constrained while spanning trees are unconstrained but global. For most practical calculations, 
like those presented here, spanning trees are more convenient. 

We want to compute one-site height probabilities, namely the four quantities P^^^h^^ = a) for 
a = 1,2,3,4. On a finite grid C, it amounts to count the recurrent configurations with a height a 
at position iq. As we will use the spanning tree description, the first task is to characterize those 
spanning trees which correspond to recurrent configurations with a height a at site iq. This was 
done for the first time in [11], and was recently reviewed, with two different derivations, in [13, 14]. 
That part is not difficult, and we only quote the results. 

Prom the burning algorithm, the spanning tree grows from its root to the open boundary sites 
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and then towards the interior of £. A site i is called a predecessor of another site io if the unique 
path on the spanning tree from i to the root passes through io (equivalently, the fire has propagated 
from io to i). Then the probabilities that the site io has height a are given by 

PciK = a) = PKK = a - 1) + (^^^' « = 1' 2, 3, 4, (2.3) 

where P^ihi^ = 0) = by definition, and where is the number of spanning trees on C where 
the site io has exactly k predecessors among its nearest neighbours (so for k = 0, io is a leaf of the 
tree). 

While the calculation of Xo is relatively easy (it was done first in [6] without the spanning tree 
description), that of the other Xi,X2 and X3 is much harder. The basic reason for this is that 
the spanning trees involved in Xq are subjected to a local constraint, while those counting for Xk, 
k > 0, must satisfy a global constraint. In the case of Xi for instance, a spanning tree is counted 
only if the firepath comes to io, then goes on straight to one neighbour of io but does not eventually 
come back to another neighbour of io after a long tour through the lattice. 

Priezzhev has devised a technique to compute the numbers X^. [11]. It turns out that they can 
all be expressed — in a complicated way — in terms of entries of the inverse toppling matrix A~^, 
that is, in terms of the discrete Green function on jC for appropriate boundary conditions. The 
thermodynamic limit of the above probabilities, when the grid C becomes infinite, then simply cor- 
responds to taking the discrete Green function on the corresponding infinite lattice. For simplicity, 
we denote the corresponding measure by P = lim^^^^ P^. 

The calculation of the 1-site probabilities, equivalently of the numbers Xk, was carried out by 
Priezzhev in the case of the infinite plane L ^1?, thereby obtaining the 1-site bulk probabilities 
Pi,P2,-P3 and P4. In the following, we carry out the calculation for the upper half-plane L 
Z+ X Z, when the boundary is either fully open or fully closed. The corresponding probabilities 
Ptt(m) depend on the distance m of the reference site io to the boundary. That dependence is 
precisely what enables us to make the comparison with field theory correlations, and to assess the 
scaling form of the random height variables 6{hig — a). 

III. HEIGHT TWO ON THE PLANE 

In this section we review Priezzhev's calculation of the higher height probabilities on the plane, 
since these calculations are the starting point for our calculations in the UHP [11]. We sketch the 
key points here; the reader is referred to Priezzhev's original paper for details. 
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FIG. 1: Labelling of sites. 



General Methods 



Majumdar and Dhar [6] had already shown in 1991 that the height one probability at a site iq 
is given by 



(3.1) 



where Xq is the number of spanning trees in which the site zq is a leaf, and M = det A (see Section 
II). 

The restrictions on Xq arc local: a spanning tree is in Xq if and only if no neighbour of iq has an 
arrow pointing to zq. Local restrictions on the spanning tree can be computed as finite-dimensional 
matrix determinants by the Majumdar-Dhar method. We briefly describe the key points of this 
method, but the reader is directed to the original paper for a detailed description [6]. 

The set of spanning trees is characterized by the toppling matrix A, which specifies the bonds 
on the lattice. When we want a subset of spanning trees with local restrictions, those restrictions 
can be imposed by changes in a finite number of entries of the toppling matrix. For example, to 
remove the bond between neighbours i and j, we simply change Aj, and Ajj to 3, and set Ajj 
and Ajj to 0. Calling the new toppling matrix A', Kirchhoff 's theorem implies that det A' is the 
number of spanning trees on the lattice with the bond ij removed, or equivalently, the number of 
spanning trees on the original lattice which do not use the bond connecting i and j. 

To compute Xq, we have to count the spanning trees in which iq is connected to only one of 
its neighbours. The simplest way to do this is to remove the three bonds connecting iq to three 
neighbours, say j2-:3d,-,3i (see Fig. 1) and to define the new matrix A' = A + B, with 

/-3 1 1 1 \ 
1-10 
10-10 

V 1 0-1/ 

where the rows and columns are labeled by io,j2,j3,j4 (all the other entries of B are zero). Then 
by symmetry, Xq = 4 det A'. 



B 



(3.2) 



10 



12 




+ 




FIG. 2: Breaking Xi into local, loop, and ^-graph terms. 



A' is an infinite-dimensional matrix, so its determinant does not make much sense. However 
the probability Pi only requires the ratio of two infinite determinants. Noting that B has finite 
rank, and defining G = A~^, one easily sees that 

det(A') 



Pi = 



= det(l + BG) 



(3.3) 



det(A) 

in fact reduces to a 4-by-4 determinant. 

It is easy to compute this finite-dimensional determinant. In the limit of an infinite planar 
lattice, G is the Green matrix of the discrete Laplacian, whose values and properties are well- 
known [28]. Using them, Majumdar and Dhar were able to show that [6] 



2(7r - 2) 
Pi = ^ ^ ' 2± 0.07363. 



(3.4) 



The other X\~ for A; > are much harder to calculate. We discuss X\ in this Section, and X2 in 
Section IV. 



B. Decomposition into local, loop, and ^-graphs 

Since all four directions are symmetric, to find we can, without loss of generality, assume 
that the single predecessor is the site to the south of zq. Fig. 1 shows the labels of the 4 neighbouring 
sites, with the predecessor site indicated by a large solid circle, and the non-predecessor sites by 
large open circles. By enumerating possible arrow configurations in Xi, Priezzhev showed that X\ 
can be broken up into the terms in Fig. 2, or the following equation: 

Xi = 12 (Mocai - ATioop + Me) , (3.5) 

where the AA's arc positive (they count certain graphs). The reader is referred to [11] for the 
derivation of (3.5), or to Section IV for derivations of similar decompositions for X2- In this 
Section, we take the decomposition in Equation (3.5) as given, and only discuss what the different 
terms mean, and how to evaluate them. 
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In the first term in Fig. 2, the arrow from ji must point to iq, and the arrow from iq must point 
to j2, while the sites and j4 are disconnected from zq. Since the conditions on the spanning tree 
are local, the size of this set of spanning trees can be calculated with the Majumdar-Dhar method 
(or Kirchhoff's theorem). The B matrix that corresponds to these conditions is 



«o 33 J4 Ji 3i-x ji+x j]_-y 



B 



/ -3 1 1 1 
1-10 
10-10 

y 0-3 







1 







1 







1 



\ 



«0 

k 



(3.6) 



Then 



local 



(3.7) 



gives the probability that a random spanning tree meets these conditions. 

In the second term in Fig. 2, the arrow from ji must point to io, and the arrow from io must 
point to _74, and there must be some chain of arrows that takes us from j4 to ji, forming a closed 
loop. The sites j2 and js are disconnected from io, and we also require that there be no other 
closed loops. Because we have a closed loop of arrows, this term is not a subset of the set of 
spanning trees. The condition on the graphs is not local, so this term cannot be calculated by the 
Majumdar-Dhar method. 

Priezzhev showed [11] how to adapt the Kirchhoff theorem to deal with closed loop diagrams 
such as these. He showed that if n bond weights in the toppling matrix are set to — e, then 

det(A') 



lim 

e— >oo 



(3.8) 



counts the total number of arrow configurations where each of the n bonds is required to be in 
a closed loop, and each closed loop contains at least one of the n bonds; in this count, each 
configuration is weighted by (—1)"^, where c is the number of closed loops. 

The second graph in Fig. 2 can thus be calculated by giving the bond from ig to ji weight — e, 
as represented by the dashed line, and making other appropriate local modifications. The matrix 
B is then given by 



12 



B 



iO J4 J3 J2 

/ -e \ 
10-10 

10 0-1 



h 
h 



(3.9) 



Since we are defining A/ioop to be the number of graphs, and each configuration gets a factor of 
(— 1)*^, A/loop/A/" is given by — det(l + BG), rather than det(l + BG), and is found to be equal to 

A/ioop _ 



1 , 2(7r-2) 

H 5 ^0,0- 



47r2 



(3.10) 



It contains a term proportional to Go,0) the inverse of the discrete Laplacian at coincident sites, 
which is infinite, see below Eq. (3.15) (for a finite lattice, Go,o diverges as logL, where L is the 
system size). The divergence of this second term in Fig. 2 reflects the fact that the number of 
diagrams with loops is much greater than the number of spanning trees. Since Xi is a finite 
fraction of the total number of spanning trees, the factors of Go,o must cancel in the second and 
third terms of Eq. (3.5). 



C. ^-graphs 

While the first and second terms in Fig. 2 can be calculated by finite matrix determinants, the 
third term is more complicated. Priezzhev called these ^-graphs. In these graphs we have two 
arrows from iQ (instead of the usual case of one arrow), to j2 and j4, and both arrows are parts of 
two separate closed loops, that go through ji to get back to io; js is disconnected from xq. The 
number Me of 0-graphs cannot be calculated by a single matrix determinant, and the bulk of the 
work in calculating P2 lies here. To do this, Priezzhev used an elegant technique, known as the 
"bridge" trick. 

The two paths from io must meet at some point before returning to io — we call that point ii. 
Priezzhev then considered two possible graphs, labelled the L and F graphs, each of which has 
three — e bonds ("bridges") between non-adjacent sites, as shown in Fig. 3. Looking at the F graph 
on the right side, we see that there are several possible ways to fulfill the condition that every — e 
bond be in a closed loop. One possibility is to make three closed loops by separate paths from a to 
J2, ii to ji, and b to ^4. The other possibility is to have one closed path that contains all three of 
the — e bonds. This second possibility can occur in two possible ways. First, by paths from a to ji, 
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FIG. 3: The bridge trick: L and T configurations. 

h to j4, and b to j2- And second, by paths from a to j^, b to ji, and ii to j2- This gives the three 
possible ^-graphs shown in Fig. 4, after reversing the direction of some of the paths, which does 
not affect the overall counting. (Note that the topology forbids cases where the three — e bonds 
only produce two closed loops. If such cases existed, they would get a different factor of (— 1)'^, 
where c is the number of closed loops, and ruin the bridge trick.) 




FIG. 4: Three possible ^-graphs from the F configuration of the bridge trick. 

When all the graphs produced by the L and F bridge configurations are considered, it is found 
that almost all the ^-graphs are listed. Characterizing the ^-graphs by the three sites adjacent to 
ii that are connected by paths to ji, j2, and j4, it is found that we are lacking all cases where the 
three sites are {ii — x,ii + y,ii — y}, and have overcounted (by a factor of 2) all cases where the 
three sites are {ii +x,ii+ y, ii — y}. However, since the ASM is symmetric under horizontal flips, 
the overcounting and undercounting exactly cancel. 
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The bridge diagrams shown in Fig. 3 are calculated with the following B matrix: 



B 



js io il a b 

/-I 1 \ 

-e 

-e 

-e 



J3 

k 



(3.11) 



Since ii can be anywhere, we need to sum over all possible locations of ii on the plane, 



(3.12) 



AfLik) and Afr{ii) are both given by det(l + BG), with S as in Eq. (3.11). Mhih) and Nvih) 
differ only in the locations of a and h with respect to i\ (see Fig. 3). While B has the same form in 
each case, the row and column indices associated with ii, a and h are different for different terms, 
which then involve different entries of the matrix G. 

However, Eq. (3.12) does not exactly give He- The bridge trick (i.e. the equivalence of the 
bridge diagram and the relevant 0-graphs) requires that the sets {iQ-,k-,k} and {ii,a,6} have no 
points in common. When ii and io arc the same, or close, Eq. (3.11) may or may not correspond 
to the subset of ^-graphs shown in Fig. 3. Priezzhev checked all such cases, and found what we 
will call "special" cases, where Eq. (3.11) either produces graphs that are not ^-graphs (in which 
case these terms need to be subtracted from the infinite sum) , or misses certain ^-graphs (in which 
cases new terms need to be added to the infinite sum). When these special terms are accounted 
for, the total number of ^-graphs is given by 



= Y. \-^L{h) +MT{ii)] +ML{k) + ^r(i4) + '^rik) -Mhih) -Mvih) 



-ML{io)-NT{io)-NT{k), 



(3.13) 



where Mrij^) is the graph obtained by adding the three — e bonds as shown in Fig. 5, and removing 
the bond between j3 and io as before. 

All of the individual terms in Eq. (3.13) can be calculated by taking finitc-dimcnsional deter- 
minants (using Eq. (3.11) for all terms except A/'t(j4))- The final difficulty is in the infinite sum 
over il. The two matrix determinants in Nhiii) and Nviii) differ only in a single column, and can 
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J3 
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• 








k 


io 


• -->- 


--• 




ii 




a 


• 







FIG. 5: The T configuration used to compute A/tUa)- The sites a and b are nearest neighbours of j4. 
thus be combined in a single matrix determinant: 



M M 



^ Gi,o — Gofi Gk,i Gk+i/ Gk,i-i — Gk,i+i ^ 

Gofi — Gifi — 1 Gk,i-i Gk,£-2 — Gk,e 

Gi^i — Gift Gk-i/ Gk/ Gk-i/-i — Gk-i/+i 

\ Gi^i — Gi,o Gk+i,e Gk+2,e Gk+i,e-i — Gk+i,e+i J 

(3.14) 



where Gf;/ stands for the entry A^j^ with {k,£) = i — j. A convenient integral representation of it 
is 



Gk 



dad/3 



_7r Stt^ 2 — COS a — cos {3 ' 
D. Final results for h = 2 on the plane 



(3.15) 



The Green function Gk,i can be computed explicitly for small values of k and I, so that the 
first column of the matrix in Eq. (3.14) contains known numbers. For the other three columns, 
the representation Eq. (3.15) can be used to write each entry as a two-dimensional integral, over 
the variables (03, /^s) in the second column, (a2, P2) in the third column, and (cti, /3i) in the fourth 
column. The sum over k and i produces two (5-functions S{ai + a2 + a3) S{Pi+P2 + P3), which allow 
to integrate over 03, /^s. So we end up with a four-dimensional integral (two further integrations 
can be carried out analytically, see Section VII). The end result is the expression obtained by 
Priezzhev [11] for the height two probability, namely 



13 2 12 1 / 4 
2 27r 7r2 7r3 4 V TT 



(3.16) 



where 



/(c„c„c3).yy_^^^yy_ 



i sin /?! • det M(ci ,02,03) 



47r2 L(ai, /3i)L(a2, P2)L{ai + 02, Pi + P2) ' 



(3.17) 
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with L(a, P) = 2 — cos a — cos /?, and 



1 



e 



1 \ 



C3 e 



M{ci,C2,Cs) = 



,j(ai+a2) 



(3.18) 



C2 e' 



1 



e 



\ C2 e 



i(ai+a2) g2ia2 




We now make an observation that will be useful for comparison with P3 later. The function 
I(,ci,C2,Cs) does not depend on C3, because all the terms in the integrand proportional to C3 are 
antisymmetric under /3i — *■ — /3i, P2 —(^2- Since the integral is linear in c\ and C2, we can write 
it as the sum of two separate four-dimensional integrals: 



Ji and J2 still need to be evaluated numerically, so this decomposition docs not help tis get a 
closed form expression for P2, but we will see in the next section that it does allow us to get a 
closed form relationship between P2 and P3. Numerical integrations yield Ji = —0.26866 and, 
quite unexpectedly, J2 = 0.5 with a error smaller than 10~^^, leading to the value P2 = 0.1739, as 
quoted in [11]. Although we could not evaluate it exactly, it is very tempting to conjecture that 
J2 is exactly equal to ^. As we will see in the next section, this would imply an exact, closed and 
simple formula for the 1-site height probabilities Pj on the plane. 



While Priezzhev calculated the height two and height three probabilities on the plane, he only 
gave the details of the height two calculation. In this Section, we give the details of our derivation 
of the height three probability, both because these details are necessary for the calculation of the 
height three probability in the UHP, and to derive an exact relationship between P2 and P3. 



For the height three probability, we need X2, the number of spanning trees where exactly two 
neighbours of iq are predecessors. This can happen in three ways, as shown in Fig. 6. Priezzhev 
showed that X^^ can be written as a linear combination of simple local graphs, as shown in Fig. 7, 
yielding (there appears to be a misprint in Eq.(32) of [11]) 



I{ci,C2) = ciJi -I-C2J2. 



(3.19) 



IV. HEIGHT THREE ON THE PLANE 



A. Height Three Decomposition 
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X. 



(1) 



.(2) 




.(3) 



O 

-o < • 



FIG. 6: Definitions of \ and \ 



X. 



(2) 
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(2) 

FIG. 7: A simple decomposition of ^ 



The calculation of X^^^ and X^^ is more complicated. Just as Priezzhev broke Xi into local, 
loop, and more complicated terms (see Fig. 2), by counting arrow configurations, we can break 
up X2'^ and X^^ into similar terms. We discuss the decomposition of Xg"*^^ in some detail. The 



derivation is shown graphically in Figs. 8, 9, and 10. 

We begin by explaining the derivation in Fig. 8. For x!^'^ , we get the same number of spanning 
trees whether the arrow from points to or j^, so without loss of generality, wc pick J4. Next, 
recall that the open circle on js indicates that it is not a predecessor of iq. We have no simple 
matrix method for calculating subsets of spanning trees with such a restriction. So we start off 
with the subset with no such restriction, and then subtract off the different ways that js can be a 
predecessor of iq: js can be a predecessor of iq either by a path through j2, or by a path through 
ji. This gives the last line in Fig. 8. 

The first graph in the last line of Fig. 8 can be calculated by the Majumdar-Dhar method. 
The last two graphs in this line require further manipulation to be put into a form that can be 



= 2 



- 2 



0-^—0 





FIG. 8: Decomposition oi X^^\ Part 1. 
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FIG. 9: Decomposition of X^ ', Part 2. 



calculated with matrix determinants. They are calculated in Figs. 9 and 10. 

In Fig. 9, we first use the fact that redirecting the io ji arrow to an is arrow will not 
affect the total number of arrow configurations. j4 is not permitted to be a predecessor of io, and 
as before, we deal with this condition by starting without this condition, and subtracting off the 
different ways that can be a predecessor of zq. This gives the second line of Fig. 9. The first two 
graphs in this line can be calculated as matrix determinants, by adding — e bonds. The last graph 
in Fig. 9 is what we call a ^-graph. In the last line of Fig. 9 we no longer draw large filled or large 
empty circles, because there are no predecessor conditions imposed, other than those that follow 
naturally from the arrows drawn. 

The analysis in Fig. 10 is similar. In the last line of that equation we again use the fact that 




FIG. 10: Decomposition of X^ ' , Part 3. 
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X. 



(3) 



= > > O- 



= 2 •-^-6-^-« = 2 



FIG. 11: Decomposition of Part 1. 



-4 






FIG. 12: Decomposition of ^ Part 2 



(3) 




dashed lines must be contained in closed loops, and the fact that directions of paths of arrows can 
be reversed without changing the total number of graphs. 

The graphical decomposition of X2 uses a similar logic, and is shown in Figs. 11 and 12. 



Defining the graphs as in Fig. 13, we have 



X. 



(1) 



8(ri-r3 + r4-r5) + i6Ar^, 



^ = 4(r2 - + 2T4) + m-g. 



(4.2) 

(4.3) 



The ^-graph is defined similarly to the 0-graph, and is shown in Fig. 13 rotated to its standard 
form. Just as with the 0-graph, there are two arrows from iq, pointing to j2 and J4, which form two 
closed loops that reach ii through ji. The only difference from the 0-graph is that now, instead 
of removing the bond from J3 to io, wc require an arrow from to i^. (The paths from J4 and j2 
now cannot go through J3, since they must still reach io via ji.) To represent the ^-graphs, wc use 
the fact that if a — e bond is added from a site to itself (i.e. along the diagonal of the B-matrix), 
it effectively turns that site into a sink, where arrows can terminate, but not begin. So by adding 
a — e bond from to itself, we fulfill the conditions necessary for a graph. 
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T 



T4 = 



y 




FIG. 13: Definitions of the Tj and ^ graphs. 



B. 0- graphs 



Since the 0-graph differs from the 0-graph only in that a — e bond has been added from j3 to 
itself, the only change from the h = 2 infinite sum is that instead of Eq. (3.11), we have 

J3 h a b 

/ -e \ 

-e 

-e 

-e 



B = 



33 

h 

k 



(4.4) 



The analysis of the ^-graphs then proceeds identically to that for the (^-graphs. There are again 
a number of special cases that arise when ii borders or is equal to io- Checking these cases by 
hand gives correction terms analogous to those in Eq. (3.13): 

^e = Y. [^L{k) + A^r(ii)] + ML{k) + Mk) + Mrik) - ^Lih) - Mh) - Ar(io). (4.5) 

n 

Nhiii) and A/'r(ii) are given by det(l + BG), where B is given by Eq. (4.4), and a and h arc located 
as shown in Fig. 3. The T graph is defined similarly to the T graph in Fig. 5, except that instead 
of forbidding the arrow from to zq, we add a — e bond from to itself. 



C. Final results on the plane 

As in the height two case, all the individual terms in Eq. (4.5) can be calculated as finite- 
dimensional matrix determinants. The infinite sum is similar to the one in Eq. (3.14). Because 
only one bond of the B matrix has been changed, only the first column of the matrix in Eq. (3.14) 
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is changed, as follows: 



/ Gifi — Go,o \ 

Gofl — Gifi — 1 
— Gi,o 



/Gi,o\ 

Go,o 
Gi,i 
\Gi,iJ 



(4.6) 



Note that this column is no longer finite; instead, all its terms diverge as Go,o- 

Just as with the height two probability, our end result for P3 is a sum of local and loop graphs, 
which can be computed exactly, and an infinite sum, which results in a four-dimensional integral. 
Since only the first column of Eq. (3.14) is changed, the four-dimensional integral still has the form 
in Eq. (3.17), and we get the same integrals Ji and J2, but with different coefficients. We obtain 
explicitly 



X. 



(1) 



AT 



= 1 - ^ + ^ - f - (^1 + ^-^2) + Go,o (4(Ji + J,) - A(. - 2)) , (4.7) 



X. 



(3) 



3 T 16 1 

= -o + r2--?-o(<^l + f-^2) + Go,0 ( 2( Ji + J2) - ^ (TT - 2) 



TT 



(4.8) 



Even for an infinite lattice, the ratios x''^^ /J\f and X^"^ /J\f must be finite since they represent 
fractions of all spanning trees. Go,o being infinite, this immediately requires that the coefficients 
of Go,o niust cancel, giving 



J1 + J2 



^(vr 



(4.9) 



While we could not evaluate Ji and J2 analytically, not even their sum, the numerical values 
mentioned before confirm Eq. (4.9). 

We find P3 from Eq. (2.3), with X2 = X^^^ + ^ + \ and then rewrite both P2 and P3 
with Eq. (4.9), reproducing the values found in [11], namely 

2 



1 1 3 12 TT 
"2 = ^-1 3- — 



2 TT TT^ TT^ 27r 



- J2 ~ 0.1739, 



1 2 12 

T + T 



8 - TT 



TT 



Ait 



-J2 ^ 0.3063. 



These two equations give us an exact relationship between P2 and P3: 

(vr - 8)P2 + 2(7r - 2)P3 = TT - 2 - - + ^ - ^. 

To our knowledge, this exact relation is new. 

Using the value of Pi = , we can rewrite the previous relation as 

48 - 127r + 57r2 - vr^ „ „ (vr - 2)(7r - 1) 

-Pi + (7r-8)P2 + 2(7r-2)P3 = ^ ^ 



(4.10) 
(4.11) 

(4.12) 



2{tt - 2) 



TT 



(4.13) 
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The significance of this new relation will be made clear in Section VIII. 

Coming back to the conjecture we made at the end of the previous section, namely that J2 = 5 
exactly (and checked to 12 decimal places), we see that, together with (4.10) and (4.11), it yields 
a very simple and exact formula for P2 and P3, similar to the formula for Pi: 

„ 1 1 3 12 „ 3 1 12 

P2 = T- ;^ 2+—' ^3 = i + 3- 4.14 

4 27r vr^ TT-' 8 vr vr-^ 

Remarkably these values imply an even simpler formula for the mean height in the bulk, 

(/i) = Pi + 2P2 + 3P3 + 4P4 = — , (4.15) 

o 

a value conjectured by Grassberger [3]. The striking simplicity of this result clearly calls for a 
better explanation than just long calculations. 

The formulae (4.14) rely on J2 being exactly equal to ^. An integral representation for J2 has 
been given at the end of Section III. By carrying out two of the four integrations, we can rewrite 

it in a simpler form, as a two-fold integral, 



J2 = — 75- 8 ^ / / sin ■ 

TT^ TT TT^ Jo V3 - COS/?i 7-7r 1 - ^1^2*3 



cos 2 COS ■ 



52 



2 



(3 — cos/3i + cos/32) COS ^ — 2sin/32 sin ^ , (4-16) 

where U = yi — ^Jyf — 1, yi = 2 — cos j3i and (3^ = —{Pi + (32)- This integral expression has been 
used for the numerical evaluation of J2, yielding J2 = 0.5 + 0(10"^^). 



V. HEIGHTS TWO AND THREE ON THE UPPER HALF-PLANE 

The calculations of the height two and three probabilities on the UHP use the same formalism 
as the infinite plane calculations. The relation between the 1-site height probabilities in the ASM, 
and the probabilities of numbers of predecessors in spanning trees, stated in Eq. (2.3), still holds 
in the presence of a boundary. The corresponding equation for the height 1, Eq. (3.1), has been 
used by Brankov et al. to obtain the 1-site height one probabilities on the UHP in the presence of 
an open or a closed boundary [7]. 

On the UHP, we will use integer coordinates (m, n), where n G Z runs along the horizontal axis, 
m G along the vertical axis, and with the boundary located at m = 1. We want to compute the 
1-site height probabilities Paim) for having a height a at the site (0, m), and when the boundary 
is either fully open or fully closed. By horizontal translation invariance, the coordinate n plays no 
role. 
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FIG. 14: Inequivalent rotations of Xi with a boundary to the south. 

When calculating Xi on the plane, we obtained a factor of 12 in Fig. 2 because of the symmetry 
between the four possible predecessors of io, and between the three possible directions of the arrow 
from iQ. In the UHP, we still have the latter factor of three, but the four choices of predecessors 
are no longer symmetric in the presence of a boundary. We thus need to break Xi up into the A, 
B, and C configurations in Fig. 14, where the dashed line indicates the boundary. (The cases with 
j2 and j4 as predecessors are still symmetric.) 

Just as in the bulk, each of these configurations can be broken up into local, loop, and 9- 
diagrams (see Fig. 2). The 0-graph can then be written as an infinite sum, and some special 
correction graphs where ii is close to io. 

The local, loop, and special correction graphs can all again be calculated as matrix determinants. 
In these terms, the only effect of the boundary is to modify the lattice Green function. The lattice 
Green functions in the presence of open or closed boundaries can be found in terms of that on the 
plane using the image method: 

^open _ ^open _^ _^ (K T\ 

^mi,m2,ni-n2 — ^(ni,mi),(n2,m2) ~ *^ni-n2,mi-m2 <J"ni-n2,mi+m2 ) W--^^ 

^mi,m2,ni-n2 — ^(ni,mi),(n2,m2) ~ ^"i-"2,mi-m2 + G'ni-n2,mi+m2-l- (5-2) 

For the local, loop, and special correction graphs, we only need the lattice Green function 
between two sites near (0, m) , and thus simply expand out the above Green functions for mi ^ 
1712 ~ m 3> 1 and ni n^- The relevant expansion in inverse powers of m reads {k.n <^ m) 

G,«„ = G„,o - ^.og(2™) - - i (I + |.og2) + '%t'l;*'^ 0(i,) . (5.3) 

where 7 = 0.57721 ... is the Euler-Mascheroni constant. 

The local graphs have the same values as on the plane, up to terms of 0{m~^). For example, 
Mocai in Fig- 2 becomes three terms corresponding to the three orientations in Fig. 14, also shown 



24 



■'^2,cl,A — 




w-loop _ 

■'^2,cl,B ~ 





FIG. 15: Rotations of graphs used to calculate Xx on the UHP. 



in Fig. 15. For the closed boundary case we have 



\ rlocal 

M 
M 



A/'2,local , 7r-2 _ TT^ - IpTT + 18 J_ 

+ , Q ,„ Q O "'"'-'^V 3/' 



47^3 



IGvr^m^ 



AA2,local ^ TT^ _ TT^ - IpTT + 22 ^ ^^ 

AAaiocai vr - 2 vr^ - Gtt + 10 



1 



(5.4) 
(5.5) 
(5.6) 



where A/2,iocai refers to the plane value in Eq. (3.7). 

The loop graphs are also matrix determinants with exactly the same B matrices as on the plane. 
Recall that on the plane, loop graphs had terms proportional to Go,0) the constant diverging part 
of the Green function. When there are boundaries, we define 

g{m) = -^[logm + 7 + ^ln2] +2Go,o = Go,o + G2m,o + C(^), (5.7) 

L ^ J TIT/ 



1 r 5 1 1 

g{m) = -— logm + 7 + -ln2 = Gq^q - G2m,o + 0{^), 



(5.8) 



to be the dimension zero terms of the Green functions in the closed and open cases respectively. 
One then finds that the loop graphs in the closed case are given by 



^loop 
Jf 



•^2,closed,C 



. rloop 

•^^ 2. closed, jj 
Jf 



1 2(7r- 2) , , 1/1 1 \ 

■4;^ + ^^^^"^^ + ^U-;?J 



+ 



1 



3 1 



13 



7r2m2Vl6 487r 2l7r2 



TT — 2 , , \ / 1 s 

g{m))+0{—), 



(5.9) 



^f 



1 I 2(7r-2) 

47^2 



g{m) 



1/11 1 



vr^m \4 



1 



1 



25 13 TT - 2 
+ :rr^ + 



7r2mHl6 487r 2l7r2 



2n 



g{m))+0{^). (5.10) 
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Once the local, loop, and special correction graphs have all been computed as series in 1/m, 
all that remains to evaluate is the infinite sum contribution to J\f{9), which is the most difficult 
part to analyze. The infinite sum contribution can be set up similarly to the bulk case, with two 
differences. First, the boundary Green functions must be used; second, the summation over ii 
must cover the UHP, rather than the whole plane. For example, for the ^-orientation of M2 in the 
closed case, we obtain the sum {{k,£) are the coordinates of ii) 



r 



/i=2,cl 



k=—oo t=l—m 



(5.11) 



with 



Ml 
N 



m+l,m+l,0 m+l,m,0 m+l,m,0 m,m,0 m+l,m,l m,m,l m+l,m,l m,m,l 



^m+i+l,m+l,k 



m-\-i,m,k 
^m+i,m,k+l 



^m+l,m,fe+l 

'~^m+i,m,k+2 
/^cl 

^m+e+l,m,k+l 



^m+e,m,k-l 
'~'m+£,m,k 
^m+i+l,m,k-l 



(5.12) 



and 



m+l,m+l,0 



*^m+l,m,0 



m+l,mfi m,m,0 m+l,m,l m,m,l m+l,m,l m,m,l 



^m+e,m+l,k 
^m+e,m+l,k+l 
^m+e-l,m+l,k 



m+£,m,k 
^m+e,m,k+l 
m+£—l,Tn,k 



^m+l,m,k+l 
^m+e,m,k+2 
^m+e-l,m,k+l 



^m+e,m,k-l 

m+£,m,k 
^m+e-l,m,k-l 



(5.13) 



The end result for the 0-graph contribution to the ^-orientation of A/2 in the closed case is 
(including the special correction graphs) 

5 2(7r - 2) , , . 1 / 1 



•^2^,closed,A _ 1 7 

M ~ 8 Svr 



471-2 ^3 

1 / 13 43 3 



-g{m) + 



3 3 \ 
2^ 16^/ 



^ m2 V247r4 487r3 327r2 ^ 647r 



m VTT 

13 (7r-l)(7r-2 



27r3 



■9{m))+lT'''''. (5.14) 



The overall result for the ^-orientation is then 

A/2.closed.A •^2°ciosed,A _|_ •^2,closed,A M'o 



,closed,j4 

77 



2, closed, A 



Af M M 

^ 1 / 1 17 17 TT - 2 N h=2,cl 



1_3_1 4 1/1_3\ 
8 87r TT^ TT^ m Vtt^ IGtt/ 

(5.15) 



Terms 4"^'^' and I'^' ' are defined similarly to ' , but with the sites in the Green functions 



r/j=2,cl _ i_r. _j „. 1 r/j=2,cl 



rotated appropriately. We then get the probability that a site a distance m from the boundary has 
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height two: 



Putting everything together, P2^{m) can be written as the sum of a power series in 1/m and 
three infinite sums of matrix determinants. It appears in the next section, along with corresponding 
expressions for other heights and boundary conditions. 

The analysis of Pz{m) in the UHP is similar to that of P2{m). The local and loop contributions, 
and the correction terms for the ^-graph, can again be calculated as matrix determinants, and 
expanded in 1/m. Just as on the plane, the key change is that for the ^-graphs, we have an — e 
bond from to itself, which changes the first row of the matrix determinants in Eqs. (5.12) and 
(5.13). 

VI. GRAPH THEORETICAL RESULTS 

We collect in this section the results obtained by putting together the various pieces correspond- 
ing to the local, loop and 0, ^-graphs contributions to the 1-site probabilities P2{m) and P3(m), 
when the reference site lies at a distance m from the boundary. For each probability, we consider 
a closed and an open boundary. As in the previous section, Ib refers to the 6-orientation, while 
Ij^C = Ia + Ic is the sum of the other two orientations. 

When the boundary is closed, the probabilities are given by the following expressions, 
1 _ TT^ - 37r2 - 47r + 24 Stt^ - 67r - 8 47r2 - 97r + 20 - 167r(7r - 2)g{m) 

- 4=^'^' - 24=2'"^ + . . . (6.1) 
pcl.^^ _ pel. X 7r3-187r^ + 96 + 247r(7r-2)ff(m) 9(7r - 4) 



47r3 167r2 



m 



. 477r^ - 1287r + 52 + 247r(20 - 77r)g(m) .=3,0! , .,/.=3,ci , 

+ ^^:;^2 +■■■ (b.2j 

with g{m) = 2Go,o - ^[logm + 7 + | log 2] = G2m,o + Go,o + 0{m-^). 

When the boundary is open, the expressions are slightly different and read 

„oDen/ X TT^ - Stt^ - 47r + 24 Stt^ - 67r - 8 llTT^ - 407r + 64 + 327r(7r - 2)d(m) 

^2^ M = ^Z3Z. + 



,op 

'AC 



- /^=2'°P - 24=2'°P + . . . (6.3) 



popen. X _ popen. x TT^ - 187r^ + 96 + 247r(7r - 2)ff(m) 3(37r^ -An- 16) 
^3 [m) - P^ [m) 

IQtt^ - 947r + 158 - 127r(20 - 77r)g(m) .=3,op ^ . r'»=3,op 



167r3 



m 



2 ^ ^^AC ^ "-'B 



+ 3/^='''°P + 6J^=^'°P + . . . (6.4) 
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with g{m) = -^[logm + 7 + | log2] = G2m,o - Go,o + 0{m-^). 

In the previous four expressions, the exphcit terms represent the local and loop contributions, 
as well as the special correction graphs encountered in the 9 and 6 graphs. All expansions are made 
to order 2 in m~^, the ellipses standing for higher order terms. Finally, the remaining Ib and Iac 
pieces correspond to the infinite sum contributions of the 9 and 9 graphs, which we now specify. 
As we have seen above, these terms include a summation over the discrete UHP, and also multiple 
integrals from the integral representation of the Green matrix entries. They are functions of the 
single variable m. 

The functions Iac in the various cases are equal to 

daid/3i e-*^"i ± gj(2m+^-l)ai rr-K dQ;2d/32 C'^"'^ ± gi(2m+€-l)a2 



lAcirn) 



If 



Stt^ 2 — cos ai — cos (3 



1 77-71 



Stt^ 2 — cos a2 — cos (32 



X // — :^^(-2zsina3 



^ Svr^ ^ 2 — cos 03 — cos /33 

where the matrices Ma and Mq take the form 



gife(^i+^2+/33) ^ip2 [det Ma - det Mc] , (6.5) 



Ma 



a+c 
2 

6+c 
2 



^ 2 



1 



gia2 
1 

gift 



1 

pi/33 



Mc 



g— c 
2 



^ 2 



3— lai 
1 

gi/3i 



1 

gi/32 



1 

gi/33 
3-2/33 y 



(6.6) 



In the integrand, the upper signs are chosen for a closed boundary, and the lower ones for an open 
boundary, independently of the height value (/i = 2 or /i = 3) one considers. The dependence on h 
is entirely contained in the values of the parameters a, 6, c, d and e, which also depend on the type 
of boundary condition: 



height 2, closed : 



o 



1 1 

'2 ^ 87rm2 ' 



c = — ■ 



2m + 1 
12 1 , ^ 



3 1 5 
+ ^ + 



4 47rm 167rm^ ' 



(6.7) 



height 2, open : < 



height 3, closed 



1 



1 

2 ~ SW 
1 



1 



2 TT Sirm? ' 



2iTm 

1 = 0, e 



1 



, , 1 1 13 

2 1 7 

6 = 2o(m) 1 \ TT. 

' TT 27rm 487rm2 



4 47rm 167rm^ ' 
2m + 1 



c = d 



47rm^ ' 



, , 1 31 
e = g(m) + h 



(6.8) 



(6.9) 



47rm 
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height 3, open : < 



-2g{m) 



1 



7 



2 487rm2 ' 



27rm 



2 1 

b = -2g{m)---——^, e = -g{m) - 
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(6.10) 



In a similar way, the other functions Ib can be written 



967rm2 



Stt^ 2 — cos ai — cos /3 



Oil /■/""' 

1 JJ-TT 



87r2 2 — cos a2 — cos /32 



//: 



dosd/^s [e*(^~-^)"3 -I- gi(2m+£-2)Q!3jg-i/32 _|„ Jg-i(^+l)a3 -|- gi(2m+£)a3jgi/32 



87r2 



2 — cos as — cos (5^ 

X e^*^(^i+^2+/33) detMs, 



(6.11) 



-i/32 



Ms 



(6.12) 



where the matrix Mb has the form 

/ a' e-*^i e-'P^ e'^'^s \ 
6' 1 1 1 

c'-d' e*"i e*"2 e*"3 

V c' + d' e-^"i e-'"-' e-*"3 / 
The choice of the signs is the same as for Iac^ and again the actual values of the parameters a', 6', c' 
and d! depend on which height variable and which boundary condition one considers, 



height 2, closed 



-- + ^ 
4 IQttw? ' 



1 



4 TT IGnm'^ 



b' 



1 



4 167rm2 



, ci' = 0. 



(6.13) 



1 1 

T - - + 



1 



height 2, open 



4 IQ-Km? ' 4 TT IQirrn? ' 

d' = 0, 



4 16irm^ 



(6.14) 



height 3, closed 



a' = c/(m) 



1 7 

+ 



47rm 967rm2 ' 



1 1 7 
- + -. + 



TT 47rm 967rm2 



I./ / ^ 1 1 1 
b' = g{m) - - + + 



height 3, open 



-g{m) 



4 47rm QGTrm^ ' 
1 



2m + 1 



(6.15) 



967rm2 



5(m) 



TT 967rm2 ' 



6' 



"^^"^^ - i + 96^^ 



1 



(6.16) 



47rm 



In order to make the expressions of P2{m) and P3{m) completely explicit, the remaining task 
is to compute the asymptotic value of the functions Iac and Ib, for the different heights and 
boundary conditions. As we will see, these functions have an expansion in inverse powers of m, 
times logarithmic corrections, which we need to compute to order m~^. 
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VII. ASYMPTOTIC ANALYSIS 



We will not give all the details of the calculation of Iac and Ib for the open and closed 
boundaries, as the analysis is technical and rather long. However, as an illustration of how this 
analysis can actually be carried out, we give here some details of the analysis of Iac for the closed 
boundary. The reader who is not interested in these technical details may skip this section and go 

straight to the next one, where the final results for the probabilities P2{m) and Pz{m) are given. 

The analysis proceeds by a long sequence of otherwise elementary steps. As a very first step, 
the similarity of Ma and Mq allows us to write the diflFerence of their determinants (appearing in 
I'^Q{m)) in a simpler form: 

c + d cos ai cos a2 cos 03 

c 1 1 1 



det Ma - det Mc 



JP2 



JI33 



i/33 



+ 



I sm ai I sm a2 i sm 0:3 
111 



JI32 



(7.1) 



c 

c e "^^ e "^^ e 

which no longer depends on e. We leave the parameters a, b, c and d free to cover the two cases 
h = 2 and h = 3. We denote by Ui and U2 the two new matrices in the previous equation, and 
accordingly we decompose Iac into two pieces, I^cim) = lf{m) + li{m). They can be handled 
in the same way, so we focus on the contribution of U2, that of Ui being simpler. 

In (6.5), the summation over k and a number of integrations can be performed: the summation 
over k produces a delta function 27r(5(/?i +/32+/?3), which then makes trivial the integration over /?3; 
the contour integrations over the three variables ai, 0:2, 03 may also be carried out. This eventually 
leaves an expression which involves one semi-infinite summation over i, and two integrations over 
Pi,/32- Before giving the result of these easy steps, we introduce the following notation. We set 



2/(/3) = 2-cos/3, m = y - - 1, u{P)=^' = y+^/f^, 



(7.2) 



and the same functions with indices yj , tj , Uj if the argument is Pj . They are even periodic functions 
of with tj (resp. Uj) being the root of — 2yjZ + 1 = inside (resp. outside) the unit circle, 
for y > 1. The range of values of y(/3) is the interval [1, 3], but more importantly for what follows, 
t{P) is strictly positive, smaller than 1, and takes the value 1 at /3 = only (or a multiple of 27r). 
The functions t{P) and u{j3) are continuous but have discontinuous derivatives at /3 = 0. 
By expanding the determinant of U2 on the first row, we obtain 



detC/2 = -ismaiD{(32,(33) + isma2D{Pi, Ps) - i sin asD {Pi, P2), 



(7.3) 
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where the function D{x, y) is a 3-by-3 subdeterminant, given by 



D{x, y) = 4i 



x-y x+y 
a cos — cos — - — 



sm ■ 



x-y 



(7.4) 



By using the following formulas, in which n G Z, 

da e*"" r da e*"" sin a 



^ 27r 2/ - cos a ^ 
da e*"" sin^ a 



1 J-n 



27r y — cos a 



zsgn(n)tl"l, (sgn(0) = 0) (7.5) 



2vr y — cos a L 2 ' J 

the summation over k and the integration over ai , a2 , and /^a yield 



(7.6) 



1 



d/3id/32 



1 



Y: [-sgn{i)tf + i^+^-i] [4^1 + ^r+^-i] [sgn(£)4^1 + 
e=i-m 



W j_ j.2m+e-li 



+ 



167r2 



1 

~ 16^ i/-. 



d/?id/3: 



V(2/i-l) (2/2-1) 

+ i(tl + i^) (t2 + t^) + + *^~')(*2 + i2"-')}, 



1 + 4""-^] 



(7.7) 



where we have kept P3 as a shorthand notation for — /3i — ^2- 

In the third double integral, which has no summation over i, the two variables can be decoupled, 
and the two simple integrals can be carried out by using 



/: 



(7.8) 



'r- 1 

The result of the integrals is then expanded for large m by using the following asymptotic expansion 
of the Green matrix {k,i <C x) 



Gx+i,k — Gx,o 



i e-k^ 
+ — ^ + ^ 



(7.9) 



2ttx 4ttx^ 

The calculations are straightforward, and one finds that the third double integral in (7.7) equals 



+ 



a - 6/1 



■ T-- + 



27rm v4 nJ IGTrm^ 



1 



/ 9 4 . ^ 



+ 



(7.10) 
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The summations in the remaining two double integrals are similar: for ^ > 0, they are equal 
up to a sign and the permutation Pi Pa; for £ < 0, they are equal up to a sign. Separating the 
three parts I > 0, £ = and i < 0, which we call respectively Yi, Y2 and Is, we can rewrite the 
first two integrals in (7.7) as 

1 /T .„ D{Pi,p2) 



+ 



167r2 



V (2/1 - l)(2/2 - 1) 

00 



2m+e-l\ _ 4.2m+e-l 



167r2 



V2/2-1 



+ 



1 



d/?id/?2 



V(y^-i)(yi-i) 



^=l-m 



(ir + *2 



t 



2m+e-l\ 



(7.11) 



the second integral in the previous expression, contains no summation, and is the simplest 
one to evaluate. If one expands the integrand, each term gives rise to a double integral which 
can be decoupled, computed exactly, and finally expanded by the procedure explained above. Yi 
and Y3 are more complicated. We will give the details for Yi only, as it contains all the necessary 
ingredients to compute not only F3 but also if, and in fact all the integrals Iac and Ib needed to 
compute the 1-site probabilites on the UHP. 

The summation in Yi is straightforward, and yields 



Yi 



1 rr 



d/3id^2 



V(2/^-l)(yi-l) 

1 (l+t2'"-l)(l+i2m-l)(^_^2m-l. 



(7.12) 



1 - tit2t3 

Because I?(/3i,/?2) is odd under Pi,P2 — A, —P2, the exponentials can be replaced by their 
imaginary parts, 

QiWi,P2) 1 



Y 



bjf 



d/3id/?2 



[1 + e^-^ + e^-^ - t"^-^ + (tli2)'"^-' - {tltzf"^-^ - (t2t3)'"*"' - {tlt2hf'^-\ (7.13) 



where the function Qi is defined as 



Qi(/3i, A) = D{l3i,l32) Jyl - 1 sin/32 - Jyl - 1 (sin/32 - sin/Jg 



(7.14) 
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In this expression for Yi , we distinguish three pieces according to the number of tj to a large power 
that they contain, namely 1, then tf^''^ + tf^'^ - 4""'^, and finally (tii2)^'""^ - (tits)^"""^ - 



A. First piece in Y^: no large power of tj 



This first term is singular: 1 / A/(yf— T)(y|— T) has a single pole on the lines /3i = and /?2 = 0, 
where Qi{Pi,f32) does not vanish (moreover (1 — tit2i3)~^ is singular at the origin /3i = P2 = 0, 
though that singularity alone is integrable). In order to isolate the singular part, we write this first 
piece as 



1 



Qi{:h:^h) Ql(0,/?2) Ql(/3l,0) 



27r2 ^(y( - _ 1) L 1 - t,t2t3 1 - tj 1 - tj 

_j_ rr d/3id/32 



(9i(0,/?2) ^gi(/3i,0) 



l-tj 



(7.15) 



Vivi - l){yi - 1) 

The second integral contains the singularity, proportional to Go,o- Using Qi(/3, 0) = Qi{0,f3) 
- ^ — 1 sin^ (3, it is equal to 

1 a-b r d/9i r sin2/32 a- 6 



d/^2 



^-(2 + -)Go,o. 



(7.16) 



The first integral is convergent by construction, and separating into parts proportional to a and 
b, can be written as —-^Ti + -^T2, where Ti,T2 are numerical integrals: 



T = r r d/32 pin(/ 

' Jo V^^-V/^^ 1- 



sin(/ ?i-/32) 
iii2i3 



y| - 1 sin/32 - \/ yf - l(sin/32 - sin/^a) 



+ 



JO \Jyi - 1 J-TT 



d/32 



- 1 sin^ /?! V2/2-I sin^ ^2 

l-tj 
sin /3i — sin P2 



(7.17) 



- 1 sin/32 - y^yf - l(sin/32 - sin/33) 



+ 



a/j/^ - 1 sin^ /3i A/yl^sin2/32 



1-t? 



+ 



1-ti 



(7.18) 



These integrals cannot be factorized, and therefore not easy to evaluate analytically, but a numerical 
integration yields the values 



Ti = 16.0897, T2 = 17.4056. 



(7.19) 



The contribution of the first term is thus 



first term of Y, = {a - 6)^Go,o - £^T, + ^r^. 



(7.20) 
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B. Second piece in Yi: one large power of a tj 

The second piece contains the three terms in Yi that involve only one tj to a large power, 
see (7.13). Making appropriate permutations of the variables and using the reflection symmetry 
-/?2) = Qi{(3i,P2), it may be written as 

(7.21) 

where QY^{x,y) = ^[Qi{x,y) + Qi{y,x)] is the symmetrized function. 

The function ti = t(/3i) equals 1 at /3i = and decreases monotonically as (3i increases towards 
TT. Therefore for m large, is exponentially decreasing, so that the dependence in m of the integral 
is controlled by the way the function / behaves near /3i = 0. 

To see this more precisely, let us first establish the following bound, which will prove extremely 
useful later, 

tiPr - e-l'^kl + ^ - ^) < + 0{m-% (7.22) 

12 96 

To prove it, one may notice that the function on the l.h.s. (i.e. the difference) vanishes at /? = 
0, and has a unique maximum which moves towards the origin and decreases in height as m 
increases. The location of the maximum, /?*, can be solved iteratively as a series in ^, yielding 
= m + + iH^ + . . .. One then simply bounds the l.h.s. by its value at /?*. The function 

g-m/3(i _j_ m| ^^) was chosen so that its Taylor expansion around zero matched that of *(/?)"* 

to order 5. 

Assume that the power expansion of f{Pi) contains a term of dimension A; > — 1, namely a term 
Pi, times possible powers of log/3i. Then at the order m~^, and for A; > —1, 

/ dAC/3?~ / d/3ie— + / d/3ie— = 0(m-i-*^). (7.23) 

JO Jo Jo 

We see that we may simply replace t^^ by e~"^^^ , and that the terms in the expansion of / of degree 
2 or higher in (3i do not contribute at order m"^. A simple dimensional analysis shows that powers 
of log Pi simply bring powers of log m in the result of the integral. 

When k = —1, the previous result essentially goes through (up to a divergent constant term 
and a logarithmic term) . Since ^y y'f — 1 = /3i + . . . near the positive origin, one has 

/•TT j.2m /•IT j.2m 

/ d/3i ^ = / d/3i J. + O(m-i) = 27rG2^,o + O(m-i). (7.24) 

Jo Pi Jo Vt/i - 1 
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Coming back to (7.21), the function /(/?i) is itself given by an integral 

^^^^^ = V-.7rT^^^^^7^FT T^pr/ = 1/^^137^- (7-25) 

To evaluate the dimension of /(/3i) near /3i = 0, one can make a Laurent expansion in /3i inside 
the integral, that is, of the integrand Gi(/?i, /32)/(l — tit2t3). One easily checks that it starts with 
a dimension —1 term near /3i = 0, and so does the expansion of /(/3i). It implies that the integral 
(7.21) will contribute to order mP (a divergent term), and to orders m"^ and m~^. 

Although this procedure of making the Laurent expansion go through the integral — we call it 
the ETI procedure. Expand Then Integrate — gives the right dimension of /, it is not valid in general 
when the actual series expansion is required. In particular, it can lead to wrong results when the 
expansion in powers of (3i of the integrand has coefficients which are non-integrable functions of 
(32, thereby producing potentially spurious divergences which the original integral does not have. 
These spurious divergences are the signal that the integral defining the function / in fact does not 
have a Laurent expansion. In our case, this phenomenon will manifest itself by the presence of 
logarithms. 

We decompose the integration domain of P2 into three pieces, [0,7r], [— /3i,0] and [— tt, — and 
first consider the region [0, tt]. The expansion of Gi in Pi is well-defined, with regular coefficients, 

Gi(A,^2) = + co{p2) + ci{(32)Pi + ... (7.26) 

Pi 

Since l — tit2t3 does not vanish (for /3i > 0), the previous expansion may be inserted in the integral 
over (32, leading to a sum of terms of the form J d(32 c(/32)/(l — ^1^2*3), where c{(32) is one of the 
coefficients in the expansion of Gi . These in turn can be written 

Th/^ _A^= r ar r c{(32) c{Q)-(3ic'{^) + (3lm/S + d'{^)/2\ ^ 
.r.o c(0)-/3ic-(0)+/?f[c(0)/8 + c--(0)/2] 

V. ^7W=i • ^ ^ 

The virtue of this decomposition is that the first integral, with its subtraction term, has a Pi 
expansion with coefficients which are integrable functions of P2 (up to order 2 in /3i), for which the 
ETI procedure may be used. Thus 

c(/?2) _ c(/32) 



1 - tit2t3 



Jo '^^ 1 — ^1*2*3 Jo 
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Prom the explicit expansion of the coefficients near the positive origin, /?2 ~ O"*", 



c_i(/?2) = -(a-6)/3i + ..., co(/32) 



a — b 



3(a - b) 



+ ..., 



(7.29) 



we obtain 



(7.30) 



The last two integrals can be computed explicitely to order in 



" d/32 

=-log/3i + |log2 + ... 



V2cos^ 

arcth- 



a/3 + cos (tt - Pi) 



(7.31) 



whereas the other one yields 27rGo,o to order 0. We find 
Jo i — 111213 Jo J- — II 



Gi(/3i^ 

1*2^3 



a — 6 

ETI 4 



/3i [27rGo,o + log /3i - | log 2] + . . . (7.32) 



The region /?2 G [— /?i,0] is left invariant under the change of variable P2 — —Pi — P2 = Pz- 
Thus the integrand may be symmetrized in P2 Ps, 



f d^2 ^^^^^"^^^ =- f d/32 ^li^liM±M^l"^3) 
J-/3i 1 — ^1*2*3 2 1 — titiH 



(7.33) 



The two variables Pi and /32 being small, a double series expansion can be used which shows that 
the integrand has dimension 4 (in the two variables). Therefore the integral over P2 brings a 
contribution at order P\ and higher, and can be neglected. 

The third and last portion [— tt, —Pi\ is sent onto [0,7r — Pi\ by the change of variable P2 ^ Pz, 
which is then further extended to [0, tt], since the integrand is of order Pi around P2 = tt: 



^1 ^ Gi{pup2) r~^'„ Giipi,p2) 



Jo 1 ~ 11*2*3 Jo 1 ~ 



*1*2*3 

where Gi{Pi, P2) = Gi{pi, Ps). 

Repeating the same calculation as with Gi, we obtain 

-01 n.(R-. R„\ r-^ h-.lR- Rr.\ 7{a — b) 



Gl{PuP2) 
*1*2*3 



(7.34) 



r\oGi{PuP2)_ r.nGiiPi 



Gl{Pup2) 
*2*3 



ETI 



/3i[27rGo,o + log/3i-|log2] + 



(7.35) 



All together, we find that the function / is given by 



f{Pi) = r dP2 f ^ - 2(a - b)Pi [27rGo,o + log /?i - | log 2] + . . . (7.36) 

Jo J- — *ll2t3 ETI 
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The remaining integral must be computed by expanding the integrand in powers of (3i, and 
then integrating term by term. Given the expansion 

Gi d_i(/32) 



1 - ^1*2^3 



Pi 



(7.37) 



we find 



Gi+Gi d_i(/32) + c?_i(-/?2) 



+ [doiP2)+doi-(32)-d'_A-P2)] 



1 - ^1^2*3 Pi 

+ Pi [dl{P2) + dii-P2) - d'oi-p2) + y-li-P2)] + 



The symmetrized coefficients read 

d_i(/32) + d-i(-/?2) = -{a-b) sin^/32 



2^2- cos/32 - Vv'i - 1 



doiP2) + doi-P2) = 



di{P2) + di{-P2) 



cos^f 



(3-cos/?2)v'y2-l ^ 



P3(cOS/32) + P2{cOsP2)Jyi - 1 



P5(C0S /32) + P4(C0S P2)\Jvl-l 



(7.38) 

(7.39) 

(7.40) 
(7.41) 



(3-cos/32)(2/i-l) 
where P2,P3, Pji and P5 are some polynomials of degree 2 up to 5. 

The integrals of these coefficients over [0, tt] can be done exactly. The first two are integrable, 
but the third one produces a divergent term due to the non-zero value of PiiX) = 4(a — h). One 
finds 

[d_i(/?2) + d-i(-/?2)] = -{a - h)^.. 



r,/^x 1, o^^ . (5^2 - 8)7r ^^2 - l)7r 
[do{P2) + do(-/?2)] = 4a + ^ - ^ 2 ' 



(7.42) 
(7.43) 



^.,ro^ r o ^^ 0/ 76 - a (39^2 - 29)7r J15V2 + 13)7r 
[di{p2) + cli(-/?2)] = 87r(a - 6)Go,o + —3— + ^ ^ - ^. (7.44) 

The integral of the derivative terms in the expansion (7.38) poses no problem and simply yields 
^^Pi- Collecting the various results, the function f{Pi) is found to be equal to 



fiPi) 



(a - 6) TT + 4 



(5V2-8)7r (V2-l)7r 



+ 4a + a- — 6^^ — + 

-log/?, + ilog2]+ 3 24 - 24 



/?i 2 ■ ■ " 2 ' 2 

3i„„oi , 76-0 , (39^2 - 29)7r ^(15V2 + 13)7r 



Pi{2{a-b) [27rGo,o 



+ a- — fe '"' } + ■ ■ ■ 



(7.45) 



It remains to multiply this function by ^ and to integrate over /3i . The integrals are carried 
out by using the following formulas: 



2m- 1 



dpi^ 



1 



Pi 



27rG2m,o + 7r- + 



7 



2m A8m? 



+ 



(7.46) 
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r - 2^ + 4^ + • • • 4^ + • • • ^'-''^ 

^ d/3ii?"*-^/3i log^i = ^[l-7-log2m] + ... (7.48) 

We eventually find 

second piece of n = \j [ d/3i t^^^VlA) 

27r^ Jo 

, ,,7r + 4^ 1 (, (5a-6)7r (26-9a)7ri 

C. Third piece in li: two and three large powers of tj 

From (7.13), the last piece reads 

J. \2m— 1 ri/'j. J. \2m— 1 (4. ± + \2m— In 



2tt^ jj_^ - i)(y2- 1) (1-^2*3) 

It is simpler to evaluate then the previous one because it contains at least two tj to a 
large power. This implies that the value of the integral is determined by the behaviour of 
the rest of the integrand around the origin, so that a double series expansion around f3i = 
^2 = can be used. The same dimensional analysis as before shows that the expansion of 



Qi'"\Pi, P2)/ \/ ivi — l)(y2 ~ ~ ^1^2*3) may be stopped to terms of global dimension (in 
the two variables). 



Special attention must be paid to the double expansions. Each of the three functions .^yj — 1 
has discontinuous derivatives on the line l3j = 0. As the integrand involves them all, it has itself 
discontinuities when one of the lines /?i = 0, /32 = and /?i + /?2 = is crossed. As a consequence, 
the partial derivatives no longer commute, and a series expansion must be computed separately in 
each of the six regions obtained when the domain [— vr, tt]^ is cut by the three lines. 

Since the integrand is even under /3i,/32 — /32, the integration can be restricted to half 
the square P2 > 0. This leaves three separate regions which we denote by (A), (B) and (C): 

(A) = [0,7r]2, (B) = < 0, /32 > 0, /33 < 0}, (C) = < 0, /32 > 0, /^a > 0}. (7.51) 

In (A), the expansion is carried out for positive values of the two variables, in any order we want. In 
(B) and (C), the expansion is made for negative values of but the order matters: the expansion 
in Pi is made first in (B), last in (C). 
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One may easily check that both the numerator Qi^™(/3i, /32) and the denominator 
\/ {y\ — l)(y2 ~ ~ ^1*2*3) have dimension 3 around the origin, in the three regions (A), (B) 
and (C). Thus the ratio has dimension zero, so that the above integral has a single contribution, 
at order m~^. The explicit expansions read 

f - &) (A + /32)(/3i - (i2? + . . . in region (A), 

Qr(A,/32) = 



\{a-h) p2{f3i - /32)(2/3i +f32) + ... in region (B), 
l{a-b) Pi - /32) + 2/?2) + . . . in region (C) , 



(7.52) 



and 



( 2/3i/32(/3i +p2) + ... in region (A), 



(7.53) 



-2/?i/?| + ... in region (B), 

^2pfP2 + --- in region (C). 

The integrals can now be computed in the three regions (the formula (7.24) is useful). The 

integration over (A) yields (the exponents 2m — 1 of the tj variables have been replaced by 2m 

since this makes no difference at order m~^) 



167r2 



d/3: 



(A) 



id/?^ ^^^^ [g-2m(/3i+/32) _ 2e-2m(2/3i+/32) _ g-4m(/3i+/32)] 



a — b 



1 5 

27rG2rrj,0 + - + - log 2 



+ ... 



(7.54) 



The other two are similar, 
a — b 

167r2 .;.;(B) 



// d/?id/32 "^^^'^^^i [e2"^(/3i -/32) - 2e-2"»^2 _ g-4m/32 



a — b 



57r^ 13 11, ^ 

- YG2,n,o + ^-^log2 



+ ... 



a — b 
"l67r2 



// 

^^(C) 



d/3id/32 



(A - /32)(/3l + 2/?2) ^^2m(/3i-/32) _ 2e2rn(2/3i+/32) _ ^4m/3i 



A/32 



a — b 



-TT^ 21 27. ^ 

647r2m2L2^2"^'° + T"T^°^^. 



+ ... 



(7.55) 



(7.56) 



Adding these three contributions, each taken twice to cover the full domain [— tt, tt]^, one even- 
tually arrives at 

a — b 



third piece of Yi 



87r2m2 



9 7 

- - - log 2 - TrG2m,0 



+ ... 



(7.57) 



D. Final results 



Adding the three contributions (7.20), (7.49) and (7.57) yields our final result for the quantity 
Yi defined in (7.13), 

y-i = ia-b)^[Go,o + G2m,o]-£^Ti + ^T2 
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+ ^^^^{7a-316-^^^^^^ + (22a-6)7r-7(a-6)log2} + ... (7.58) 

As mentioned earlier, the other terms Y2 and ^3 can be computed in a similar way, with 
analogous results. Together with (7.10), Yi, Y2 and ^3 complete the calculation for if- The 
evaluation of if proceeds the same way and, with if, yield I^q- The three other integrals, namely 
-^B' ^AC ^^'^ ^ handled by the same techniques. The calculations become fairly long, but 
each step is straightforward if one follows the methods used in the sample calculation detailed 
above. The next section summarizes the final results. 



VIII. SINGLE HEIGHT PROBABILITIES ON THE DISCRETE UHP 

In Section VI, explicit formulas were given for the single site probabilities P2{m),Pz{m) on the 
UHP, at a site located a distance m from the boundary, chosen to be open or closed. We have 
shown in the previous section how to evaluate the asymptotic expansion of the four functions, I'^Xci 
-^S' ^'ac -^s'' second order in m~^. We now put the various pieces together and give the 
final results for the probabilities. The divergences found in the individual terms, like in (7.58), all 
cancel out, as they should, leaving finite probabilities. 

We obtain for the closed boundary, 

1 / 34 , 4(7r-2)/ , 5,_ N , 4(7r - 2) 



P|'(m) = P2-7r^<^-ll + — + ^ ^ 7+-log2 +^ ^logm^ (8.1) 



Svr^m^ 1 TT TT V 2 



X „ 1 f 5 44 2(8 - 7r)/ 5 , \ 2(8 - tt) , 1 

and for the open boundary, 

= + ^ ( -9 + - + (7 + I log 2) + log ™1 , (8.3) 



St^'^w? \ tt tt V 2 

„oD. X „ If 3 36 2(8 - tt)/ 5, \ 2(8-7r), 1 

= ^^+8^r+2-v+^^(^+2i°^^)+^^^°^4- (^-^^ 

The first terms, independent of m, are the bulk probabilities, and so can be written in terms of 
the results from Section IV. 

Prom the above expressions, it is not difficult to see that the probabilities take the form 

Pf{m) = Pi - (oi + 6i log m) + . . . , (8.5) 
ifP(m) = Pi + -^{ai + I + 6, logm) + . . . , (8.6) 
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with the following exact values for the coefficients aj and bi, 

7r-2/ 5, \ ll7r-34 , tt - 2 



£12 



27r3 



/ tl , \ liTT - 34 ,77 — 2 , „, 



as 



» - / 5 , „\ 271^ + 57r - 88 ^ 8 - tt 
-4;^(^+2^"g^)+ 16.3 ' ^3 = ^. (8.8) 

For completeness, we recall that the 1-site probabilities Pf{m) and P°^(m) have the same form, 
with the following coefficients [7], 

ai = bi = 0, (8.9) 

while those for P^^{m) and P^^{m) are obtained by subtraction, 

04 = -(«! + 02 + 03), b4 = -ibi + b2 + b3). (8.10) 

As we shall see in the next sections, the identity ai = 62 is not accidental, but follows from the field 
assignment for the height 1 and 2 variables, namely the height 2 scaling field is the logarithmic 
partner of the height 1 field. Similar relations will hold for 63 and 64. 
The numerical values of these constants are 

ai = 0.0184091, 02 = 0.0402789, 03 = -0.0154395, 04 = -0.0432485, (8.11) 
h = 0, 62 = 0.0184091, 63 = 0.0391728, 64 = -0.0575818. (8.12) 

The values of 02 and 62 have been announced in [24], though not in an exact form for 02. In the 
same reference, numerical values had also been given for 03,63,04 and 64, and obtained by fitting 
numerical simulations of the sandpile model with an Ansatz for the height 3 and 4 conformal fields. 
The values quoted in [24], namely 03 = -0.01243, 63 = 0.03810, 04 = -0.04636 and 64 = -0.05667 
(04 and 64 were obtained from an independent fit, not by subtraction), are in good agreement with 
the exact values given above. 

The forms of (8.5) and (8.6) make it manifest that the probabilities Pj(m) — Pi for i = 3,4 are 
linear combinations of those for z = 1,2, at order m~^, since they are all linear combinations of 
the two independent functions and m"^ log m, 

Psim) -P2 = as [P2(m) - P2] + A? [Piim) - Pi], (8.13) 
P4{m) -Pi = a4 [P2{m) - P2] + p4 [Pi{m) - Pi], (8.14) 

for the two boundary conditions. The coefficients a, and /3j satisfy = aia2 + PiUi and bi = aib2 
for i = 3, 4, and take the actual values 
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TT -I- 4 32-1- 4-Tr -I- TT^ — TT*^ 

a, = -3.12791, /34 = 4.49453. (8.16) 

Remarkably, it follows from these values that the linear relation we have derived above in (4.13) 
for the probabilities at infinity, actually holds everywhere in the UHP (for large m and at order 

m~^, that is, in the scaling limit) 

48 - 127r + Stt^ - TT^ „ ^ ^ ^ ^ , , „ / x (7r-2)(7r-l) 

Pi(m) + (tt - 8) P2(m) + 2(7r - 2) Psim) = ^ ^ ^. (8.17) 

z(^7r — ■i) TT 

The constant part, corresponding to m = oo, is the relation we had derived at the end of Section 
IV. It is truly remarkable that the coefficients of the non-constant parts, proportional to 1/m^ and 
logm/m?, obtained after a long and painful asymptotic analysis, so conspire to fit the same linear 
relation. 

In fact, we argue in the next section that the relations (8.13) and (8.14) actually hold at the 
level of the fields representing the four height variables, so wc expect the last formula to hold in all 
generality, that is, in any geometry and for any sort of boundary condition. Multisite probabilities 
would be linearly related in the same way, leaving the probabilities for heights 1 and 2 as the only 
independent ones (again in the scaling limit). The physical origin of this relation remains to be 
understood. 



IX. CONFORMAL FIELDS 



The conformal field theory interpretation of the previous results has been given in [24], which 
resulted in the field identification of the height variables in the scaling limit. We recall here the 
salient features of this identification, give computational details that were omitted in [24], and 
make further comments. 

The conformal field theory believed to describe the scaling limit of the sandpile model has 
central charge c = —2. It is non-unitary, and logarithmic [25], meaning that the theory possesses 
reducible but indecomposable representations of the Virasoro algebra. This characteristic property 
implies that correlation functions contain logarithms. 

We first review the basic ingredients we need to analyze the lattice results given in the previous 
section, referring to [26] for a more complete description of the c = —2 conformal theory (see also 
[29] and [30] for more general reviews on logarithmic conformal theories). 

The 1-site probabilities Pi{z) = {5(h{z) — i)) are expectation values of the random variables 
6{h{z) — i), computed with respect to the stationary measure of the sandpile model. In the 
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continuum limit, the subtracted random variables S{h{z) — i) — Pi, normalized to have a zero 
expectation value on the plane, are expected to converge to conformal fields hi{z, z), such that the 
multipoint field correlators reproduce the scaling limit of the multisite probabilities. 

The lattice results show that the four fields hi have a scaling dimension equal to 2. The height 1 
field hi{z, z) was the first one to have been studied [6-9], and turns out to be a primary scalar field 
(p{z,z) with conformal weights (1,1). This field is a true primary field, with rational multipoint 
correlators. 

It has been argued in [24] that the other three fields h2{z,z),hs{z,z) and hi{z,z) are all asso- 
ciated with another field ^(z, z), of conformal dimensions (1,1), the logarithmic partner of (f){z, z). 
The fields (f) and -0 form a rank 2 Jordan block under dilations since their transformations read 

Lo^ = (p, LoV = V'-i0- (9-1) 
The precise statement made in [24] was that the lattice calculations for the height 2, as well as 
the results of the simulations for the heights 3 and 4, now evidenced by the exact results presented 
above, are all compatible with the following assignments, 

hi{z,z) = (p{z,z), h2{z,z) = il){z,z), (9.2) 
hi{z,z) = azil){z,z) + I3i(t){z,z), h4{z, z) = a4il;{z, z) + P4(f){z, z), (9.3) 

with the coefficients ai,(3i given in the previous section. 

The general form of the 1-point functions of the fields hi{z, z*) on the UHP match the dominant 
terms of the 1-site probabilities Pi{m) — Pi, when z — z* = 2im. The normalizations of hi and 
/i2 can be uniquely fixed by comparing with, for instance, the closed boundary results for Pi and 
P2. However we have to make sure that the same normalizations would be obtained had we chosen 
to make the comparisons with the open boundary results. This is a highly non-trivial consistency 
check because much more of the structure of the c = — 2 theory is needed to prove that indeed the 
same normalizations would follow. 

The probabilities PP^(m) and P°^{m) refer to a different boundary condition, and correspond- 
ingly the two 1-point functions {hi{z, z*))'^' and {hi{z, z*))°p are related to each other by the change 
of boundary condition, itself implemented by the insertion, in the correlator, of a specific bound- 
ary field /X [31]. In the present case, we know from [20] that the field which effects a change of 
boundary condition from closed to open, and vice-versa, is a primary chiral field with conformal 
weight h = —1/8. 

It may therefore be used to compute the exact relationship between the open and closed results, 
and verify that it agrees with the lattice results (8.5) and (8.6). The general form should be retained. 
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but the coefficients of l/w? and logm/w? should be related as (— aj,— (oj + ^,6i)°P. This 
check probes correlators with two fields ji, and therefore the operator product (OPE) of /x with 
itself. This OPE contains another logarithmic field w, which itself proves extremely important to 
make the conformal picture fully consistent with the lattice results. 

An even stronger and independent support for the above field identifications will be provided by 
comparing conformal predictions with actual numerical simulations of the sandpile model. Using 
again the jjl field, one may compute the expectation values {hi{z, z*)) on the UHP with the negative 
real axis closed and the positive real axis open. They may then be conformally mapped onto an 
infinite strip with open and closed boundary conditions on either side of the strip, and eventually 
compared with the 1-site probability profiles across the strip, measured from numerical simulations 
on a finite but sufficiently long strip. This is done in the next section, where we will observe a very 
good agreement between the two sets of curves. 

In order to relate the open and closed boundary expectation values of the fields /ij, we consider 
the more general situation of an open boundary condition on the real axis with a finite stretch 
/ = [-2;i,2;2] which is closed. The corresponding 1-site probabilities P°'^'^{z), which are no longer 
translation invariant in the horizontal direction, ought to be given in the scaling limit by 

- F. = = '"'^''y'f '■'fr"»°' . (9.4) 

{H{Zi) H{Z2))op 

The notation (. . .)op implies that the boundary condition far out, namely around the point at 
infinity, is open. Hence the product i-i{zi) n{z2) really means /Lt°P'*^'(2;i) fi'^^''^^ {Z2) , since the closed 
stretch is in between zi and Z2 ■ The all open or all closed cases is recovered upon taking the limit 
zi, Z2 ^ ov —zi, Z2 +00 respectively. While the denominator is simple and can be taken to be 
(/Lt(zi) fJ,{z2))op = 2^12^5 the numerator is more complicated and is computed in the next subsection. 



A. Conformal correlation functions 



The non-chiral 3-point function is equal to a 4-point chiral correlator, in which two chiral fields 
hi are inserted at z and at the mirror image point z*, namely the left chiral field at z, and the right 
chiral field at z* [31]. As hi is a combination of (p and tp, we have to compute the chiral 4-point 
functions {^ficpcj)) and {^iiipip), the latter requiring in any case the former. So we concentrate on 
{lx{zi)ii{z2)ip{z2,)ip{zA)). 

Because the field jJL is degenerate at level 2, the chiral 4-point function can be computed in 
the usual way, by solving first the conformal Ward identities, and then an ordinary second order 
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differential equation. A first complication arises because the field not only is the logarithmic 
partner of a primary field (f), but is also not quasi-primary. Under special conformal transformations 
(Li and Zi), it transforms into the fields p and p, of conformal weights (0,1) and (1,0) respectively. 
We will take the OPE of t/' with the left and right stress-energy tensors to be 

T{z)tp{w,w} = -3+ "2 + + ... (9.5) 

[z — wy [z — w)'^ z — w 

T{z)ip{w,w) = -ri rT3+ jz + — — + ■... (9.6) 

[z — wy {z — Wj'^ z — w 

The factor —1/2 in front of (f) is conventional (when if) is normalized to be the height 2 variable in 
the sandpile model, it ensures that (f) has the right normalization to be the height 1 variable). For 
more generality, we replace it by A in what follows. 

A second complication comes from the fact that the fields p and p are themselves not fully 
primary: p is a left primary (with left weight 0), but is not quasi-primary under antiholomorphic 
conformal maps, 

T(.W„,^) = + (9.7) 

f ») = + + M^:^ + . . . (9.8) 

[z — wy [z — wy z — w 

It is the converse situation for p, which is right primary but not left quasi-primary, 

, , , K p(w,w) dp(w,w) , , 

T{z)p{w, w) = ^ + ' + ^ + . . . 9.9 

[z — wy [z — wy z — w 

r^, ^ , N dp(W,W) 

T{z)p{w,w) = +••• (9.10) 

The coefficients of the cubic terms are equal (to k), on account of Lip = LiLiip = LiLitp = Lip. 

The inhomogeneous transformations of V' implies that the correlator {ij,fj,ip{z)tp{z*)) depends on 
eight other 4-point functions, in which at least one tp has been replaced by either ^ or p. (Note 
that the fields iIj{z) resp. ip{z*) are subjected to left resp. right conformal transformations only, 
so that the chiral p and p behave as primary fields; since p behaves under right transformations 
as p does under left transformations, we use the same notation p for both.) As p and cf) are chiral 
primary fields, the correlators {fifipp), {iip,4)p), (fifipcf)) and {fifKpcj)) can be computed separately 
as solutions of a second or first order homogeneous differential equation. Once these are known, 
the four correlators containing one ip field can be computed, namely {pppil^), (ppipp), {ppip4') and 
(ppcpip), which depend on the first four. Finally (ppipip) can be computed as a solution of a second 
order inhomogeneous differential equation. More details on the procedure may be found in [29]. 
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As usual, the conformal Ward identites are used the fix the kinematical factors in the correlators, 
leaving an unknown function of the cross ratio x, which we will take to be x = The Ward 

identities read ((...) is any one of the nine correlators) 

{J2di}{...)=0, (9.11) 

i 

[J2{zidi + hi) + di]{...) = 0, (9.12) 

i 

{ ^(z,^a, + 2h^Zi) + 2zik + . .) = 0. (9.13) 

i 

The operator 5 and e take into account the inhomogeneous parts in the conformal transformations 

of V' under Lq and Li: 5ip = X(p and iip = p, whereras 6 and e annihilate all the other fields. 

All nine correlators contain a fi field and as a consequence of its degeneracy at level 2, they 
satisfy a second order partial differential equation. 



1=2 

A correlator which contains a (f) satisfies a different second order differential equation, since (f) is 
also degenerate at level 2. If the field (j) is at Za, it reads 

Idl _ y \^L_ + -h±%- + -^^1 )(...) = 0. (9.15) 

As they may be of interest for other calculations, we give some details about the derivation of 
the relevant 4-point correlators. The easiest correlators are those without V'- As their calculation 
is standard, we merely give the results, 

(/x(l)/x(2)p(3)p(4)) = Fpp{x), = Ai + A2 log (9.16) 

1 + V 1 ~ X 

1/4 _ 2 

(m(1)M2)</>(3)0(4)) = ^F^^ix), F<^^ = 5^=, (9.17) 

^/'^ h _ 
(m(1)m(2)<^(3)p(4)) = F^,{x\ F^, = (9.18) 

^13 ^23 X 
^1'^ / I _ 

(M(l)Ai(2)p(3)0(4)) = ^^F,^{xl F,^ = D^L_^, (9.I9) 

■2-14 -224 X 



The five coefficients Ai,v42,S,C and D are integration constants. 

The conformal Ward identities satisfied by the correlator (n^ijjp) require that it be of the 
following general form, 

5/4 

(/x(l)/i(2)V'(3)p(4)) = \f^, - XF^, log (^i^) + ^ fJ . (9.20) 

^13^^23 L ^ Z12 ^ Z12 -I 
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When one inserts this form into (9.14), one finds that the unknown function F^p satisfies an 

inhomogeneous ordinary differential equation, 

7x 1 — oc 1 — oc 

x(l - x)Fl + (3 - + -^F^p = 2(1 - x)F'pp + -^Fpp + XF^^, (9.21) 

where we have used the equations satisfied by Fpp and F(f,p to simplify the right-hand side. The 
general solution involves two additional arbitrary coefficients, 

-1 2\C^/T^. 



F^p{x) = E,^^^ + E2{ 



2 

- + 

X X 



X. 1-^/r 

-log- 



+ 



2 

X 3 a; 



1 + yr^ 



1 ^ 1- 

logx- - log 

^ 1 + — X 



X 



■ \ogx 



— x^ 

— .T J 



(9.22) 



The 4-point function {fifipip) is related to the previous one by the permutation 3 <-> 4, which 
induces on x the transformation x 



. One thus obtains 



(/.(l)/.(2)p(3)V(4)) 



5/4 
Z14Z24 



- XFp^ log ( — )+ — Fpp 

\ Zi2 / 212 



(9.23) 



and 



Fp^{x) = Fi 



+ F2 



2(1 -x) VT^, 1-VT^-i 2XD^/T^ 1-x 



+ 



log 



+ 



+ 



A2{ 



log- 



X 



20:^, 1-x 1 1-vT^ 
+ 3^^"S^-2^°S + vT^ 



2 I. X 

The next case is the 4-point function {^ij,ip4>). Prom the Ward identities, its general form is 



(9.24) 



(/i(l)/i(2)V'(3)</)(4)) = 



"12 



'34 



/'Zi3Z23\ 



F^^ - XF^^ log j 



X Z34 



P0 



(9.25) 



2:12 / X — 1 224 

When it is inserted into (9.14) and (9.15), one finds that satisfies two second order ordinary 
differential equations. By combining them, one obtains a first order equation, namely 



2(1 - x)(2 - x)F;^ - xF^^ = -lil^±^^^AF^^ + ^(4 - 5x + x')Fp^. 



The general solution reads 
Ff4,{x) = G 
Similarly, one finds 



x-2 UXB-D) . 2XB + D x-2 , 

, V 1 — X , logx. 



(/x(l)Ai(2)0(3)V'(4)) 



1/4 

h2 



'34 



F^ - XF^^ log J + X— F^p 



223 



with 



^ , , X-2 UXB + C) 2XB -C x-2 1-x 
F^{x) = - \ ' + ^= log ■ 



Vl — X 3^/l — X 



(9.26) 



(9.27) 



(9.28) 



(9.29) 
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We finally come to the last correlator {fi/jtipij;) . Prom the Ansatz 



(/x(l)/x(2)V'(3)V'(4)) = 



1/4 

h2 



'34 



+X^Fss Uzi) + XF^p Mzi) + XFp^ Mzi) + Fpp fs{zi)\ , (9.30) 



we find that the Ward identities are satisfied for the following choice of the functions fji 

/l(..)=log(^), /2(..)=l0g(^), /3(..) = ^^, = -X^, (9.31) 

V -212 / V Z^o ) X-lZoA Zl?. 



k{Zi) = h{Zi)f2{Zi), h{Zi) = fl{Zi)fi{Zi), friZi) = f2{Zi)f3{Zi), h{zi) = 

Inserting this form in (9.14) yields an ordinary differential equation for F^^, 

33/ , I X 



Z23 
'34 



-213-^14 



(9.32) 



x(l - x)F;^ + (1 - — )4^ 



X 



2(1 -x)2 



F.,pp + 



2(1 -x)^^" 
2(4 - 4x + 5x2; 



3x(2 - x)2 



2 - .T 
-X'F^ 



-XFa, 



(1 -x)(2-x) 



, 2(4 -3x) 2(4-x) x(6-6x + x^ 2x2(2-x) 



2(1 - x)2 1-X 



3(l-x)(2-x) '"^^^ 3(2 -x) 

in terms of the eight functions derived earlier. 

The general solution to this equation would be fairly complicated to write down and is in any 

case not needed. We are primarily interested in the above correlators when the two /i fields form a 
pair n^t',^- 

(zi)/x'=^'°P(z2). When zi goes to Z2, the operator product expansion closes on fields which 
live on an open boundary. As shown in [22], it reads 



M°P'"1(zi)m^1'°P(z2) = Z,^^"I + 



1/4 - 



(9.34) 



Therefore the limit zi2 cannot produce any logarithmic singularity, whatever the other fields 
in the correlators are. Imposing this condition on the first eight correlators, brings five constraints 
on the coefficients, namely 

X^B 



A2 = 0, C = -D = -XB, E2 = -F2 



(9.35) 



These relations reduce the number of arbitrary constants in the right-hand side of (9.33) from 11 
to 6, namely Ai, B, Ei, Fi, G and H. The most general function F^^ contains the general solution 
of its homogeneous equation, and also six functions which solve the inhomogeneous equation, each 
one corresponding to one of the six coefficients remaining in the r.h.s. of (9.33), 



Fipip — F^ijj 



homog 



+ F, 



~t~ F^^ -\- F^^ 

Ai B 



El 



+ F, 



tptp 



+ F, 



+ F, 



H 



(9.36) 
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These seven functions are given explicitly by 



r X-2 
homog y i — X 

_ 4 - 4a; - 



x-2 , 



:log 



4(1 - x) ' 

l={l + (.-2)logi^}, F, 



x/r^ 1 x-2 



Si 3vT^ 



Fi 



logx, 



=^|2-3x + 2(x-2)log^|, 

1 — X L X J 



G 30~- X - X 

|x + 2 - 2(x - 2) logxj, 
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(9.41) 



2 - 16 log 2 2 - 2x + x^ 4 X - 

9 3(1 -x) 9"7r 



X 



{Li.(- 



= ) +Ll2(- 



1 + Vl 



— X J 



2, 1 — X 10 r logx , , X ~| 8 X — 2 , , x 



Ir^ ^-2, /l-xl-0:^\-i ^ 1 



9^1^ 
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(9.42) 



where the polylogarithmic function Li2{z) = p- is regular at the origin. 

Imposing now that the correlator {fifi^tp) itself is free of logarithmic singularity at Zi2 = 0, and 
moreover that liuiz^^^o Zi2^^ (/x(l)/x(2)'i/;(3)^(4)) depends on 2:34 only, leads to five new conditions 
on the coefficients, 



B = 0, Ei = -Fi 



-Ai, H = G, l2 = -(^i-AG), 



(9.43) 



leaving only three free parameters, Ai,G and Ii. The explicit solution is 



{H{zi)n{z2)tp{z3)tp{z4)) = 



1/4 
h2 



"34 



{( 



^ Ai 2XG\ x-2 ^i4-4x-3x2 
^ ~ ~3 3~/ ^ 1 1 -X 



--(Ai-AG)^_log 
3 vl — a; 



(l + VT^)^ 4, 



X 



_ (^,_AG)-2AG.^ 

V 1 — o V 1 — X 



^ ( ZO/I \/l — 



X 



/ 2 Vzo.'i 



X 



22:34 \ 



Z2W^ - X Z23 J''"" ' 2 \Z23 {x - 1)Z24^ ZiaZu)^^^ 



log 1^34 1 



}. (9.44) 



We have required that the correlator be regular when 21,2:2 — ^ 0, namely when the boundary 
is all open. The other extreme regime corresponds to an all closed boundary, obtained by letting 
—zi,Z2 +00. In this limit, the correlator should reproduce the 1-site probabilities in presence 
of a closed boundary and thus should also be regular. This is the case provided 



A^ = AG. 



(9.45) 



At this stage, the correlator {ij.ij,ip{z, z*)) = {iJ.iJ,ip{z)'tp{z*)) physically relevant to our situation 
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depends on the two arbitrary coefficients Ai and Ii, and reads explicitly, 

{n{zi)fi{z2Mz3,Z4)) = ih - Ai) + 2Ai log [2:341 

2:34 -^1 — X 4 L — X -x/l — X 

\Z2W^-X Z23 ' 2 \Z2Z (X-1)2;24 2^132:14^ 

Anticipating a little bit, we trade the coefficients A\ and I\ for two other coefficients and 62, 

^1 = -62, h = 2a2 - b2{^ + 2 log 2). (9.47) 



In terms of those, the correlator becomes 

.1/4 

-J 260 Inp- ■ ^ 
1 1 



x-2 r 

(/x(2:i)/x(22)V'(2;3,2:4)) = /- 2^2 log 

I4 



+ 62 



^34 



62 62 — 2 
+ 2a2 + -^- 



+ 



2 4 

(9.48) 



' -213-^24 \-X-2 2V1 - X- 

The other correlator wc need involves the non-chiral field z*) = j LqiP{z, z*) = j LqiIj{z, z*), 
which becomes (j){z, z*) = (j){z)ijj{z*) = 'ip{z)(p{z*) in terms of chiral fields. The relevant correlator 
is thus given by (9.25) or (9.28) with B = C = &nd H = G = ^ = namely 

(M(Z1)/X(22)0(Z3,Z4)) = ^4=- (9-49) 

^ ^34 s/i- — X 

When we will come to the comparison of the conformal predictions with the numeri- 
cal simulations, wc will also need two other correlators, namely {fi{zi)ii{z2)p{z3, Z4^)) and 
(/x(zi);u(2;2)^(2;3, 2:4)). In terms of a pair of chiral fields, p{z,z*) = Litp(z,z*) is identified with 
p{z)'tp{z*). The corresponding correlator is given in (9.23) and reads (we use the values of the 
constants found above, and Ai = —62) 

{piz,Mz2)p{zs,Z,)) = -b2zli'{ + + V-]- (9-50) 

Z\iZ2A X ^Z\i ^Z2A ^ 

Similarly, identifying the field p{z,z*) = Liip{z,z*) with ^(2)^(2*), the correlator in (9.20) yields 



{p{zMz2)p{z3,Z4)) = -b2zli'{ - + T- + T-}- (9.51) 

^13^23 X 2213 2223-' 

The value of k in the transformations of p and p can also be fixed. The field k1 = Lip{z,z*) 
should be associated with p{z)p{z*), so that the correlator (9.16) yields 

K{p{zi)p{z2)) = -b2zli\ (9.52) 
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B. Closed versus open UHP probabilities 

Let us now come back to the 1-site probabilities on the UHP to see how they are related when 
the boundary is either fully open or fully closed. We set A = — ^ in the expressions of the previous 
subsection. 

We first consider the height 1 probabilities, already discussed in [22]. According to (9.4), the 
identification hi{z,z*) = (p{z,z*) and (9.49), the probability Pi^'\z) is equal to (2:4 = = z*) 

P^^'\z) - Pi = z;,'/' (m(zi) m(^2) Hz, z*)) = j^^^ (9.53) 

When zi and Z2 go to 0, the cross ratio x = ^V'Zlyi^Z^l goes to 0, while y'l — x goes to +1. 
In this limit, P'^'^ {£) — P\ must go to the open boundary result. This fixes 62 = = '5'' 

lim P°P-^(z) -Pi = P°P(m) -Pi = (z-z* = 2im). (9.54) 

If now —z\^Z2 go to +00, the cross ratio still goes to but \J\ — x goes to —1 because \ — x 
circles around zero when —z\ = zi take the values from to +00. Therefore we obtain the expected 
change of sign, 

lim Pl^^Hz) - Pi = Pfim) - Pi = — ^. (9.55) 

-Zl,22— >+00 

The probability P'^'^ {z) that the site at z has height 2 in presence of an open boundary except 
on the interval / = [2:1, 2^2] which is closed, is equal to the correlator (9.48) divided by z!^^ (there 
was a misprint in the last term of the formula (8) of [24]) 

z - z* 



{z-z-Y I ^ ^ 2 '^'24 



+ 2a2 + — - 



X 



{z\- z)[z2- z*)\.x -2 2Vl-a::-'-' 
with X = (^^I^)(^^Z^*) as before (although a; is a complex number, x* = the probability is 
real). 

Upon taking zi and 2:2 to 0, the expression reduces to 

lim P°P'^ (z) - P2 = ^ [a2 + ^ + 62 log m] = P°p (m) - P2 , (9.57) 

which confirms that the two coefficients 02,62 are those derived in the lattice calculations, and 
given in (8.7). It is now an easy check to see that the other limit -2:1, 2:2 +00, bearing in mind 
the monodromy factor of \/l — x, yields exactly the result for the closed boundary, namely 

Um P^P'-^ (z) - P2 = - [02 + 62 log m] = P^ (m) - P2 . (9.58) 

-2:1,22^+00 ^ 
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The same limits also hold for the height 3 and 4 probabihties. The relations (8.13) and (8.14) 
expressing P3(m) and P4{m) as linear combinations of Pi{m) and P2{m) suggest that the fields 
h^i^z) and h^^z) are given by the same hnear combinations of hi{z) and h2{z), as in (9.3). Then 
the formula (9.56) holds for P°^'^(z), i = 1,2,3,4, provided the coefficient a2)^2 are replaced by 
ai,bi, given in Section VIII. Consequently the open and closed probabilities P°^{m) and P^^{m), 
for z = 3 and i = 4, are correctly related to each other by the insertions of /j, fields. 

So the lattice results for the open and closed boundary conditions are exactly compatible with 
the effect of a change of boundary condition computed in the conformal setting. The formula (9.56) 
is however far more general, and gives access to a probability that looks at present very hard to 
obtain from lattice calculations. 

Beyond the two extreme cases of an all open or all closed boundary, this formula may be 
confronted with numerical simulations in an intermediate situation. Before doing this, we return 
to the pure closed boundary with a further comment. 

Relying on the operator product expansion (9.34), we have imposed that the correlator 
(jLt(l)/Lt(2)'0(3)V'(4)) should not produce logarithmic singularities in the limit zi,Z2 — 0. On phys- 
ical grounds, we have then required the same regular behaviour in the limit —zi,Z2 — ^ +oo since 
the 4-point correlator reduces to the closed boundary probabilities. In that limit, the two fields ^ 
come close to each other at infinity, so that the operator product expansion may again be invoked. 
However the relevant expansion is not the same as in the other limit, because it is now the open 
stretch which is reduced to nothing, in such a way that the expansion must close on fields which 
live on a closed boundary. The proper operator product expansion reads [22] 

fi^''°P{zi)ii°P'^\z2) = zli" [lon{z2) - - I log Z12] + ... (9.59) 

TT 

where the boundary field lon has conformal weight 0, and forms with the identity, a rank two 
Jordan block under dilations. The interpretation of lo{z) in the sandpile model is that it introduces 
dissipation at z [22]. 

The insertion of the previous expansion in the 4-point function yields in the limit 

(V(z,z*))ei= lim \{u;^(z2mz)^iz*))-^^^{^P{z)^P{z*))}. (9.60) 

The absence of logarithmic singularity then implies that the chiral 2-point function {iIj{z)iIj{z*)), 
appropriate for the closed boundary, vanishes identically, but also that 



{iIj{z,z*))c\ = {un{oo)iI){z)iI){z*)). 



(9.61) 
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Though strange at first sight, this relation has a natural and physical origin. Probabilities on 
infinite lattices are limits of probabilities calculated on finite lattices. The sandpile model is not 
well-defined if the lattice does not have dissipative sites. In the model considered here, dissipation 
takes place at the open boundary sites, namely at those sites which arc connected to the sink, 
itself located all around the lattice. In the case of the UHP, three boundaries have been sent 
off to infinity. When the remaining boundary (the real axis) is open, the sandpile model is well- 
defined without introducing further dissipation on the other three boundaries. When it is closed, 
dissipation has to be introduced on the other three boundaries by connecting them to the sink. 
The sink is, with the boundaries, subsequently sent to infinity, which is where the dissipation is 
eventually located. This is exactly what the field un{oo) does: it inserts at infinity the dissipation 
needed to make the sandpile model well-defined. 

A similar relation holds for the (p field, describing the height 1 variable, 

{^{z, z*))a = {cjNi^)(l>{z)cl>{z*)). (9.62) 

The insertion of ujn prevents {(j){z, z*))ci ~ Pf{m) from being identically zero, since the chiral 
2-point function of (p vanishes. 

X. FINITE SIZE CORRECTIONS VERSUS SIMULATIONS 

The simplest case after the pure open and pure closed boundary condition is when the boundary 
is half closed and half open. It is also a case which is relevant for a comparison with numerical 
simulations. 

We take the boundary of the UHP to be closed on the negative side of the real axis, and open 
on the positive side. It clearly corresponds to taking zi = — oo and Z2 = 0. Prom (9.56), we obtain 
for i = 1,2,3,4, 

^*'<^' - - -w^ {-< - l^h - M 

The conformal map w = ^ log z transforms the upper half plane to an infinite strip of width 
L. The negative and positive real axis are mapped respectively to the line = L and = 0, 
which are therefore respectively closed and open. 

In order to transform the probabilities from the UHP to the strip, we need the conformal 
transformation law of the logarithmic field (the field (p is primary). The integration of the 
infinitesimal conformal transformations (9.5) and (9.6) yields 

1pstrip{w,w) = \z\w)f^1p^i,p{z,z)-^log\z'{w)\'^4>^i,p{z,z) 
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+ 



z"(w) 



z"(w) 

Puhp{z, Z) + „ Pnhp{z, Z)+ K 



z"{w) 



2z''^{w) 



;io.2) 



where the last piece comes from the transformations^ of p and p. 

Taking the expectation value, one obtains the corresponding 1-site probabilities on the infinite 
strip in terms of those on the UHP and the expectation values of p, p. For the height 1 probabihty, 
the relation is 



(10.3) 



whereas for the height 2 probabihty, 



Pf '^iw) - P2 



z'{w)\'[[Pt\'>'{z) - P2] -'-\og\z'{w)\' [Pf °P(.) - Pi] 



z"{w) 



z"{w) 



/cl|op 



}. (10.4) 



2z''^{w) 

The height 3 and 4 probabilities are simply obtained by linear combinations (see (9.3)). In these 
formulas, the expectation values refer to the UHP with the boundary condition closed on the 
negative axis and open on the positive axis, i.e. z\ — — oo and Z2 — 0. Prom (9.50) and (9.51), 
these are given by 



(p(^,^*))ci|op = ^2 { 



^/zfi 



z — z^ 



+ 



{p{z,z*))a\ 



°P "^\z-z 



* ^ 2z/' 



(10.5) 



and (I)ci|op = 1. 

Putting everything together and setting u; = u + w, we find the following 1-site probabilities 
on the strip, for z = 1, 2, 3, 4, as functions of v only. 



jstrip 



TT ^2 COs{tTv/L) 

sm^{7rv/L) 



(i) 



r .bi 



1 + cos(— ) 

1j 



+ 



Mo.(^sM^))}.|r. (10.0) 



5,,7r2 



Should the field identifications proposed earlier be correct, these formulas represent the scaling 
limit of the probability that a site of an infinite strip, located at a distance v from the open 
boundary and L — v from a closed boundary, has a height equal to i. 

We have run numerical simulations on strips of width N and height M for various values of M 
and N (details about these simulations have been given in [24]). For ratios ^ large enough, we 
expect that the results on a finite height strip be very close to those on an infinite strip, so that 



^ The inhoniogcncous terms in the non-chiral transformation laws of p and p have been overlooked in [24] , resulting 
in a number of missing terms in the equations (10-12) and (14) of that reference. In particular the equation (14) 
of [24] misses a constant term proportional to l/L^, see the equation (10.6) in the present article. This however 
makes no numerical difference. 
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FIG. 16: Subtracted probabilities to find a height 1 (orange), 2 (in blue), 3 (green) or 4 (red) at a site lying 
along the short medial line of a rectangle M x N , and going from an open boundary (left end of the graph) 
to a closed boundary (right end). The dots correspond to data obtained from numerical simulations, with 
M — 200 and N = 50, while the solid curves represent the conformal predictions (10.6), with L = + i, 
for the limit M = +oo. 

the numerical simulations should be reliable enough to be compared with the exact scaling limit 
(10.6), provided by the conformal approach. For simplicity, we chose the boundary conditions on 
the two sides separated by M to be open — exactly what boundary conditions are imposed there 
should make no difference — , and the other two were open and closed. The two sets of plots, the 
numerical results and the conformal formulas, in each case for the four heights, are shown in Fig. 
16 (as explained in [24], the width L in the continuum and the width N of the lattice should be 
related by L = N+^). The simulation data are the same as in the figures of [24], but the conformal 
curves are plotted here with the exact coefficients Oj, 6j, derived from the lattice calculations, rather 
than with fitted coefficients. 

The agreement is very satisfactory for the four heights. In the previous section, we had shown 
that our proposal for the field identifications were fully compatible with the lattice results for 
homogeneous boundary conditions. Here we have an independent and direct confirmation. 

XI. NATURE OF THE LOGARITHMIC FIELD rp 

The previous sections contain two major results: first, the computation of the 1-site probabilities 
for having a height 2, 3 or 4 on the UHP, and second, the identification of the scaling limit of the 
corresponding height variables with a logarithmic field ■0) belonging to a logarithmic conformal 
field theory with c = — 2. There exists a well-known free field realization of a local logarithmic 
theory with c = —2, so that a natural question is whether our field ■i/' can be realized in this free 
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field theory. Our answer to that question will be negative. 



A. Comparison of correlations 

The free field theory defined by the following action 




(11.1) 



where 6, 9 are anticommuting scalar fields, is a logarithmic conformal field theory, with central 
charge c = —2, and is a realization of the so-called triplet theory [25-27, 32, 33]. Zero (constant) 
modes, if any, play a special role, due to their grassmannian nature. Because the action does not 
depend on them, their integration "eats up" one pair of 0, 6 fields in the expectation values if there 
is such a pair, and yields identically zero if not. Apart from this peculiarity, correlators of local 
functionals of B and B can be computed in the usual way from Wick's theorem. 

Correlation functions are usually normalized by the partition function, with the result that 
(I) = 1. When the fields 9,9 have zero modes, the partition function vanishes identically, so that 
this normalization is not possible. A convenient alternative choice is to set {ui{z,z)) = 1, where 
u = uq -.99 : + (3q \s a. dimension zero scalar field, which forms with the identity I a rank two Jordan 
cell, 

Loa; = -^L (11.2) 

On the other hand, when there are no zero modes, the usual normalization (I) = 1 can be chosen. 

In [22], the bulk field u and the boundary field to^ describing the insertion of isolated dissipation 
at a bulk site and at a site of a closed boundary respectively, have expectation value equal to 1, 
{ijj{z, z))piane = {'^n{x))c\ = 1- They can be realized in terms of the free fermions as 

a;=-l:^^:+-l(^+^log2) + l, a;iv = +- log2. (11.3) 

Their transformations under dilations are Lqu = — and Lqun = — ^I- 

On the infinite plane, the fields 9, 9 have zero modes, ^ and ^. The above normalization implies 

(I) = 0, {9iz, z)9{w, w)) = m = 27T. (11.4) 

Higher correlators (and OPEs) are computed using the following elementary Wick contraction, 

9{z,z)9{w,w) = — log\z — w\. (11-5) 
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For instance, the 4-point function on the plane takes the form [27] 

{9{zi)e{z2)9{z^)e{z^)) = 27r[log \zi - Z2\ + log \Z2, - Z4\ - log \zi - Zi\ - log \Z2 - Z2\]. (11.6) 

On the UHP with closed (Neumann) boundary condition on the real axis, namely 89 — BO = 
89 — 89 = on M, the fields also have zero modes. Thus the relations (11.4) still hold, while the 
Wick contraction appropriate to the closed boundary condition is 

9{z,z)9{w,w) = -log\{z - w){z - w)\. (closed on M) (11-7) 

If the open (Dirichlet) boundary condition is imposed, ^ = ^ = on R, the fields no longer have 
zero modes. In this case, one has (I)op = 1 and the 2-point function {9{z)9{w))op is simply given 
by a Wick contraction, equal to 

9{z,z)9{w,w) = -log -. (open on M) (11.8) 

I I Z UJ 

The same is true for the mixed boundary condition we have considered earlier, namely closed 
on IR_ and open on IR+. The fields do not have zero modes, and the relevant Wick contraction 
reads (the phases of ^/z, \pw are between and '/r/2) 

9{z, z) 9iw,w) = - log (v^- Vw){V^+V^) _ ^^j^gg^ ^pg^ ^^^ 9^ 

I 1 (v^- vw)(v5 + v^) 

Let us now come back to the 1-point functions on the UHP we have computed earlier for the 
height variables, to see what they are in the free fermionic theory. 

The scaling field cf) corresponding to the height 1 variable has been much studied. It is a primary 
field with conformal weights (1,1), realized in the free theory as [8] 

(j){z,z) = -Pi -.8989 + 8989: = -Pidd :99: . (11.10) 

In order to compute the expectation value of (p on the UHP, one may use the OPE to make a point 
splitting, 

:99: = lim \9{z)9iw) +\og\z - w\], (H-H) 

and then evaluate Wick contractions. 

For an open or a mixed boundary condition, the 1-point function of 4> is equal to 

(</>(z, z*))op = {4>{z, z*))ei|op = -Pi dd hm [9{z)9{w) + log \z - w\]. (11.12) 

For the closed boundary condition, the boundary dissipation field u}n{oo) is to be inserted, but is 
eaten up by the path integration over the zero modes (the piece in a; at proportional to the identity 
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does not contribute since 4> contains derivative fields only). Once the zero modes are integrated 
out, one simply compute Wick contractions for what remains, so that the previous formula is also 
valid for the closed boundary. 

Using Eqs. (11.7)-(11.9), one easily obtains 



{w(oo)«2,z-))d = ^j^^, «2.^')>op = -(j^^, (11.13) 



«.,.-))*p = -fp|i^. (11.14) 

in agreement with our earlier results, Eqs. (9.54), (9.55) and (10.1) with a\ = Pi/4 and h\ = 0. 

The scaling fields corresponding to the heights 2, 3 and 4 involve the logarithmic partner of the 
primary field (p. In the fermionic theory, the logarithmic partner of (p is proportional to 

Mz,z) = : ee{dede + Bede) ■.+idd:ee:=:eeddee:+idd:ee:, (11.15) 

and may contain an arbitrary multiple / of its primary partner. The expectation value of this term 
is given by the expressions we have just established. 

The same point splitting procedure as for ^ can be applied to ipio, for instance, 

^o(z) = - lim [9d9{z) 9de{w) + ede{z) ede{w) - singular terms]. (11.16) 



The subtraction terms are such that the 1-point function {tp0{z, z*)) is equal to the limit w ^ z of 
the point splitted fields {9d9{z) 9d9{w) + 9d9{z) 9d9(w)), in which one removes the singular part 
in each Wick contraction. 

Again for an open boundary or a mixed one, the 1-point function follows from Wick contracting 
the point-splitted fields, resulting in 

{i^eiz, ^*))op = (^e(-2, ^*))ci|op = ^™ { [e(^{w) + log \z - w\\ d9{z)d9 {w) 

- 9(z)d9(w) - -—^ — -1 \d9{z)9(w) + — A + cel. (11.17) 

For the closed boundary, the insertion of a;jv(oo) = ^ :99: +^Ilog2 brings two contributions. 
One comes from the identity piece, and reduces to the calculation we made above for the field 
(f). If we remember that the integration over the zero modes equals 27r, one obtains that this first 
contribution is 

2 41oB-2 

-log2(V.,(.,z*))ei = -^^^. (11.18) 

The second contribution from the :99: piece produces logarithmically divergent terms, because 
after the integration of the two zero modes, one 9 (or 9) in (jJn{(X)) may be contracted with a 9 (or 
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9) in ^5t. On the other hand a regularized, finite answer is obtained if we decide to integrate the 
two zero modes contained in the :66: piece of ujn{oo), since this leaves us with contractions of 9, 9 
within ^pg. This regularized expectation value is then given by the same formula as for the open 
or mixed boundary, Eq. (11.17), and yields, by using the appropriate Wick contraction, 

(^^^)cl,reg = {i + log |Z - (11.19) 

Computing the required Wick contractions for the open and mixed boundaries, and adding the 
contribution from the term I dd :99:, we obtain the 1-point function of tpe on the UHP for the 
three boundary conditions, 

{Mz, z*))op = (^^-^ + log 1^ - ^1 }' (11-20) 

(wjv(oo)V'e(^, ^*))cl,reg = ^ { - / - 4 log 2 + i + log |z - | } , (11.21) 

We see that the first two expectation values of V'e not related as they are for the field 
associated with the height variables h > 2. In particular, they do not show the characteristic 
change of sign of the logarithmic term. Moreover, the third one {''Peiz, z*))c\\op does not have the 
form obtained in (10.1) for {ip{z, 2;*))ci|op- This suggests that the fields ijj and ijjg are distinct fields, 
and that ip has no realization in the free fermionic theory. 

A further check of this statement is as follows. The two fields ijj and ipo are different, but 
they have the same abstract conformal transformations. Namely, the OPEs with the stress-energy 
tensors (9.5) and (9.6) are formally the same for ip and ipg. Consequently, the correlation func- 
tions {iJ,{zi)fi{z2)ip{z, z*)) and {^{zi)ii{z2)^e{z, z*)) satisfy the same differential equations, so that 
{lJ,{zi)iJ,{z2)ip0{z, z*)) must be of the same general form (9.44), as derived in Section IX.A. 

The general solution (9.44) involved three arbitrary parameters, Ai,G and Ii. For the tp 
corresponding to the heights 2, 3 and 4, we had required that the limit —zi,Z2 — ^ oo, in which 
the boundary becomes fully closed, does not produce the logarithmic singularity {'tp{z, z*)) log zu 
coming from the identity channel in the OPE ij,{zi) iJ,{z2) , 

hm ^-//'(/x°P'-'(2i)/x^i'°P(22)V'(^,^*))= lim^ {{uNiz2Mz,z*))-^-^^{iP{z,z*))}. 

(11.23) 

This had implied {ip{z, z*)) = 0, and the relation Ai = AG. This condition cannot be imposed for 
4'e because the expectation value {ipe{z, z*)) is not zero on the UHP. 
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Thus for il)Q, the limit in (11.23) does produce a logarithmic singularity, but this singularity is 
in any case discarded since, according to our prescription of inserting a field ujn{oo) for the closed 
boundary, we are only interested in the ujn channel of the OPE, namely the first term in (11.23). 
As it turns out, the limit (11.23) also produces the divergent terms logjzi — z\ and log|z2 ~ -2*1) 
which correspond to the divergencies discussed above in terms of the free fcrmions. The rest, which 
forms the regularized expectation value, should be a function oi z — z* , which, in order to conform 
with the results for {un {00)1^0 {z, z*))ci and {tpe{z, z*))op established above, ought to contain a term 
log I2; — 2;*| with the same sign as for the open boundary (corresponding to the limit Zi,Z2 — ^ 0). 
It is not difficult to see that this condition imposes the relation Ai = —2 AG, reducing the general 
solution (9.44) to 

1/4 



,,,,,,, , , z,i) ( ^ X - 2 AiA-Ax - . x -2 

(M.i)M^2)^.(.3,.4)) = ^{^i7r^ + T^3^i ^^7r^^°s 



(l + vT 



y/l - X 4 1 - X - X 234-^/1 - X 

+ Ai 7^-- -^34 + ^ 7 r^ + ^^W34-2Ai . 11.24 

\Z2/^\/\-X Z2Z / 2 VZ23 {X-1)Z24 Zi^ZuJ J 

The limit zi, 2:2 — > yields, after the division by .2;i2^, the expectation value of tpo in front of an 
open boundary, 

{Mz, z*))op = f ~\.2 - log 2) + + ^1 log \z - z*f\ . (11.25) 

(Z Z ) I J 

The limit —zi,Z2 +00, after discarding the singular terms log j^^^ and dividing by ^12'^, yields 
the regularized expectation value in front of a closed boundary, 

(Mz, ^*))ci,reg = , \ - 2(/i - 2Ai log 2) - 8^1 log 2 + Ai + Ai log \z - z*\^] . (11.26) 

{z — z*y L J 

Both expressions coincide with the results (11.20) and (11.21) obtained in the fermionic free theory 
if we choose Ai = —\ and I = —2{I\ — 2Ai log 2). 

Moreover, the limit zi ^ — oo and Z2 —>■ yields the expectation value of i/^e in front of a 
half-closed half-open boundary. 



*^2 /ii^i2 



/ , / *NN Z + Z^ , 4:Z{Z-Z^) ^ lZ-Z^r-4:\Zr ^ „ . , 1 

\z\{z-z*Y^ (V^-I-Vz*) 4\z\{z + z*) ) 

(11.27) 

which, remarkably, exactly matches the form (11.22) obtained in the 9, 9 theory, for the same values 
of Ai and /. Again we conclude that the logarithmic field ipg in the fermionic theory and the field 
ijj describing the height variables h> 1 correspond to two distinct fields. 
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B. Representation theory 

In the previous subsection, we have shown that the field if) describing the heights /i > 2 in the 
sandpile model has no realization in terms of free fermions. On the other hand, the fields ^ and 
uj corresponding to the height 1 variable and to the insertion of dissipation in the bulk, do have 
such a realization. Likewise the fields p and p appearing in the conformal transformations of ip can 
also be realized in the free fermionic theory. So all the fields we need except the ^ip can be realized 
in terms of 9, 0. This peculiar situation can be understood from a representation theory point of 
view. Indeed the conformal transformations typical of these fields allow for an infinite number of 
inequivalent representations [34]. We show here that the two fields -0 and il)s in fact belong to 
inequivalent representations. 

Let us first examine the situation in the free fermionic theory. For the sake of clarity, we attach a 
suffix 9 to the fields constructed out from 9 and 9. So we take ipe = 99{d9d9+d9d9) and its primary 
partner (j)Q = d9d9 + d9d9. The Wick contraction of 9 and 9, 9{z, z) 9{w, w) = — log — i/;!, implies 
the normalization of the stress-energy tensor, Tg = 2 8989. 

Prom the OPE of tpe with T, we obtain 

Loi^e = V'e - ^(pe , Liipe = pe, (11.28) 

with pg = ^{989 - 989) a field of weight (0, 1). There is a second field pg = \{989 - 989) = Li^g 
arising from the right conformal transformations of ipg. All these fields belong to a reducible but 
indecomposable representation TZ, which is however a representation of an extended algebra [26] . 
Clearly the two fields pg and (j)g are not independent but satisfy 



L-ipg = 




(11.29) 



The keypoint of this section is to show that, although a relation between pg and (pg is to be expected 
in general, the value of the coefficient ^ is particular to the triplet theory, of which the symplectic 
fermion theory is a realization. 

In Section IX, we have computed the correlator {p,{zi)p{z2)il^{zs, Z4)) where V' is a field with 
the same conformal transformation as ipg. The general form of this correlator has been obtained 
on the sole basis of the conformal transformations of the fields involved, namely (the first relation 
had been written LqiP = %p + X(p but A was taken to be equal to — | later on) 

LoV' = tp-^(f>, Liip = p, Liip = p, (11.30) 
Lqp = 0, Lip = 0, Lop = p, Lip = kI, (11.31) 



61 



and similar relations for p (moreover L_i and Z_i act by differentiations). The important fact in 
these calculations is that we have not used the identification of L_ip with a multiple of ^. 

The general solution for {ii{zi)ii{z2)i){z'i, Z4)) is given in (9.44). It depends on three ar- 
bitrary coefficients, Ai, G, /i. The derived three-point correlators {^{zi)^{z2)p{z^, Z/^)) and 
{fi{zi)^{z2)(l){zz-, z^)) are given in (9.23) and (9.28) respectively, where we can use the values of 
the coefficients as given in (9.35) and (9.43). One obtains 



5/4 rr— . 

(m(-^i)/x(z2)p(z3,Z4)) = ^i^l^i ^ + ^ + ^1 — |, (11.32) 

Z14Z24 ^ X Z Z\2 J 

1/4 



x-2 



{l^{zi)f^{z2)cf){z3,Z4)) = G (11.33) 



To see whether can be consistently identified with ^, we simply compute the partial 

derivative of the first correlator to obtain 



4i z^l^ T-2 

as (m(^i)//(^2)p(^3,^4)) = ^ ^^5=- (11-34) 



A simple comparison with (11.33) implies that the general solution allows the following identifica- 
tion 

L_ip = ^</. with /3 = ^- (11-35) 

In the interpretation of the height 2,3 and 4 variables of the sandpilc model in terms of the 
field, we have set A\ = AG = — ^ in order to avoid a logarithmic singularity when the two \i fields 
are sent off to infinity (see (9.45)). This gives 

/3asm = -\. (11.36) 

In the previous subsection, where we made contact with the fermionic theory, we argued that 
the correct choice to reproduce the fermion correlators is A\ = —2 AG = G. This then gives 

Pe = \, (11-37) 

in agreement with the above relation (11.29) obtained from straight calculations. 

Since the normalizations of p and are completely fixed by that of V', the parameter (3 cannot be 
scaled away, and this implies that different values of (5 label inequivalent conformal representations 
[34]. So there is a one-parameter family of inequivalent Virasoro representations from which a 
field V' transforming as above can be chosen. The specific if) that is relevant to the sandpile model 
belongs to the representation with /3 = — i, while the one realized in the free fermion theory 
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corresponds to /3 = 5. Because the coefficient G can be traded for /?, the general solution of 
{^{zi)fi{z2)'tpi3{z3, Z4)) for a fixed given value of /3, depends on two arbitrary coefficients Ai,Ii, as 
expected (one for the norm of ip, the other for the multiple of 4> that can be freely added to it). 

This explains why our general correlators involving p, p and (j) (and lo) arc separately consistent 
with what these fields would be in the fermionic theory. But when wc bring in the it auto- 
matically adjusts the relative normalization of p and in a way that is not consistent with the 
fermionic picture, except if the parameter (3 is fine-tuned, and taken equal to ^. 



In this last Section, we build on the knowledge we have of the scaling description to assess 
further properties of the lattice height variables in the sandpile model. We will focus on the 2-site 
correlations for height variables on the plane, and deduce their large distance limits. On the lattice, 
these correlations have never been computed, so that the formulas which follow can be considered 
as predictions, or conjectures. 

If Pij{zi2) denotes the joint probability that the height at zi be i and the height at z-2 be j, 
the 2-site correlations are given by Pij{z\2) — PiPj for i,j = 1,2,3,4. In the scaling regime, when 
the distance zu is large, these correlations should be equal to expectation values of pairs of fields 
hi{z\, zi)hj{z2,Z2)- However the plane has no boundary where sand can leave the system, so that 
the prescription we used earlier requires to insert a bulk dissipation field a;(oo) at infinity. Thus 
one should have, in the scaling limit, that 



As the four height variables hi{z, z) are linear combinations of the fields (/) and ip, we clearly need 
to compute the (non-chiral) 3-point functions {ipipio), {(pipto) and {(pcpto). We concentrate on the 
first one, and let first the uj field be inserted at Z3. 

The general form of {iI>{zi,zi)iIj{z2, Z2)io{z3,zs)) can be computed from the Ward identities, but 
the nested procedure as we used it in Section IX, is rather cumbersome. Instead it is simpler to use 
the finite Mobius transformations directly, by seeing the insertion points zi,Z2, zs as the images of 
a fixed triplet, say wi = —1,W2 = +1, = 0. 

The conformal tranformation u> — *■ z of the plane onto itself which brings {wi,W2,W3) = 
(— 1,+1,0) onto three arbitrary points izi,Z2,Z3) is 



XII. BULK CORRELATIONS OF HEIGHT VARIABLES 



Pij{zi2) - PiPj = {hi{zi,zi)hj{z2,Z2)uj{oo)), 



Z12\ » 1- 



(12.1) 



w{z) = — 



Zl2{z - Zs) 



(12.2) 



Z23iz - Zi) + Zisiz - Z2)' 
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under which the conformal transformations of tp and lo read 



\w' {z)\'^^'il^{w, w) — - log \w'{z) 



[w, w ) 



w"{z) , , w"{z) , , 

+ IZ-^TT^ P{W, W) + _ p{w, W)+K 



2w;'2(z) 



2m'^{z) 



w"iz) 



2w'^{z)\ 



"}■ 



u}{z,z) = u){w,w) — —\og\w' {z)'\^ . 



(12.3) 
(12.4) 



It is straightforward to compute the various derivatives, and express ^(2:1,^1) as a hn- 
ear combination of ^(— 1, — 1), 0(— 1, — 1), p(— 1, — 1), /J(— 1, — 1) and the identity, and similarly 
for ip{z2,Z2) and a; (23,^3). Inserting these expressions into the 3-point function yields 
{iIj{zi, zi)ip{z2, Z2)lo{z3, zs)) in terms of fifty constants, given by various correlators of fields in- 
serted at —1,-1-1,0. In the limit 2:3 — 00 however, most of these fifty constants combine or drop 
out. Specifically, when Z3 ^ 00, the transformations of the fields read 

4 



ip{zi,zi) 

1p{z2,Z2) 
^(^3,^3) 



{V(-l) + log <^(-l) - P(-1) - P(-1) - 

.(l)+p(l)+p(l)+Ac}, 



V'(l) + log 

1 



^(O)-^log 



Z12 



(12.5) 
(12.6) 
(12.7) 



'2.Z13Z23 

Using these, one sees that the 3-point function retains a logarithmic singularity in the limit 
2:3 ^ 00, unless the following identities hold (among others) 



{'lp{zi,Zi)tl){z2,Z2)) = {(t){zi,Zi)'tl){z2,Z2)) = {(f){zi, Zi)(f){z2, Z2)) = 0. 



(12.8) 



The 3-point function then takes the form 

16 



{tl){zi, Zl)'ll}{z2, Z2)UJ{00)) 



Z12 



+ log </'(-l) - P{-1) - P{-1) - K 



1) + loe 



Z12 



{!) + p{l) + p{l) + n oj{0) 



+ 2B log 1^12! + ^log2|zi2||, 
where the coefficients A and B are given by 

A = 16(<^(-l)<^(l)a;(0)), 

B = 16((/.(-l)[V(l) -0(l)log2 + p(l)+p(l)]a;(O)). 
The same procedure with one ^ replaced by a <p, or the two ijj's replaced by ^'s, yields 



(12.9) 
(12.10) 

(12.11) 
(12.12) 



{(p{zi, Zl)jp{z2, Z2)uj(oo)) 
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V(l) + log I ^ I </>(!) + P(l) + P(l)] cu(O) 



= T^\B + Alog\zi2\}, 



(12.13) 
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and 

16 A 
{cp{zuzi)cl>{z2, Z2M^)) = I ^ (</)(- l)</>(l)a;(0)) = , ^. (12.14) 

From these results and the way the height variables are expressed in terms of (p and ijj, namely 
hi{z, z) = aiilj{z, z) + (3i4>{z, z), one easily obtains the scaling form of the 2-site correlations, 

Pijizu) - PiPj = -^{aiajC + {ai(3j + Piaj)B + PiPjA 

\Zl2\ ^ 

+ [2aiajB + {oiPj + (3iaj)A] log \zi2\ + aiajAlog^ |2i2|}- (12.15) 

The value of A is related, through (12.14), to the correlation of the two height 1 variables on 
the plane, whose dominant term was computed in [6] , 

^iiW = ^'-§ + - (12-16) 

This fixes the constant A = — Pf/2. The other two constants B and C cannot be determined 
purely from conformal theoretic arguments, and require further input from lattice calculations. 
The value of A alone however determines the dominant terms of all correlations, and yields 

Piiir) - PiPi logr, p,.[r) - PiPj ^ -^^^^ \og\, i,j = 2,3,4, (12.17) 

where the values of aj have been given in Section VIII, namely q;2 = 1, as = 2.12791 and a4 = 
—3.12791. The heights are found to be anticorrelated except the height 4 variable, which has a 
positive correlation with the other three (but is anticorrelated with itself). 

In order to check these behaviours and to assess the values of B and C, we have run further 
numerical simulations of the model. The main problem here is that the correlation term ~ 
in the joint probabilities Pij{r) is largely subdominant with respect to the trivial constant piece 
PiPj- It means that we cannot test these scaling behaviours in the regime where they are supposed 
to hold (r S> 1), since that would require a huge numerical precision on the data. We have 
nevertheless tested the scaling formulas for small and intermediate values of r (r < 20) , and found 
very good agreement, in some cases for as small distances as r = 2. Also the fitted values of B, C 
are consistent for all fits. 

We have run the simulations on a cylindrical lattice, of height L = 50 and perimeter P = 200. 
We have sampled a total of 3.5 x 10^ recurrent configurations, and for each configuration, we have 
measured the heights on the six circles located at mid- height. On each circle, and for each value of 
r between 2 and 20, we have recorded all pairs of heights which are r sites apart, and subsequently 
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computed the occurrence probabilities Pij{r). Thus r is a distance which was always directed along 
the perimeter of the cylinder. 

Once the joint probabilities Pij{r) are obtained, we fit them (leaving out the data for r = 2, 3) to 
the form given by (12.15), that is, a constant contribution PiPj and terms log'^ r/r'^, for /c = 0, 1, 2 
depending on the values of Then, using the exact values of the coefficients Ui and (5i, the 
values of PiPj,A,B and C could be extracted (note that, as a further check, we did not use the 
known value of PiPj and ^4). As expected, better results are obtained for larger values of The 
results for P13, P33, P34 and P44, as well as the fitted curves, are shown in Fig. 17. 




5 7.5 10 12.5 15 17.5 20 5 7.5 10 12.5 15 17.5 20 



FIG. 17: Two-site correlations of height variables on a cylindrical lattice (height L = 50, perimeter P = 
200), measured along the perimetrical circles located at mid-height. The abscissa represents the separation 
distance, in lattice sites. The dots correspond to data obtained from numerical simulations, and the solid 
curves represent the fitted functions, given in (12.15). 

As a first qualitative check, we see that the signs of the correlations are as predicted, namely 
negative, except for the height 4 with a strictly smaller height (P34 here). 

As mentioned above, the agreement is excellent, and is so, somewhat unexpectedly, down to 
very small values of r (r = 4 and even smaller). For Pi^^r), the fitted values of the subtraction 
term P1P3 and of A and B were found to be (no log^ r term was included in this case) 

P1P3 = 0.02261, A = -0.002648, B = -0.004861. (12.18) 

The first two are consistent with the exact values P1P3 = 0.02255 and Agxact = —0.002711. 

exact 

The fits of i-*33(r), i-34(r), P44(r) yield the following values, 

P3P3 = 0.09410, A = -0.002821, B = -0.004718, C = -0.008883, (12.19) 
P3P4 = 0.1366, A = -0.002517, B = -0.004329, C = -0.008773, (12.20) 
P4-P4 = 0.1982, A = -0.002699, B = -0.004065, C = -0.009608. (12.21) 
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Again the values of the subtraction terms compare well with the exact values, P3 



0.09381, P3P4 



exact 



= 0.1367, Pi 



exax:t 



exact 



= 0.1991, and the fitted values for A are also compara- 



ble to the exact value quoted above. The values of B and C obtained from these four (or three) 
independent fits are all close to each other, from which we infer that the exact values should be 
around 



B ~ -0.0045, C ~ -0.009. 



(12.22) 



As far as we can see, these values can only be confirmed by exact calculations of correlations on 
the lattice. 



XIII. CONCLUSIONS 



Let us first summarize our main results. As a preliminary but crucial step, we have carried out 
the calculation of the 1-point probabilities of the four height variables on the discrete upper-half 
plane. This was long and technical, but gave us an exact formula for their scaling term 
where m is the distance to the boundary. As 1-point functions on the upper-half plane are 2-point 
correlators, this enabled us to obtain our second main result, namely, the field identification of the 
four bulk height variables in the scaling limit, in terms of conformal fields of a c = — 2 logarithmic 
conformal theory. The picture that emerges is that the height 1 variable turns out to be a primary 
field, whereas the heights 2, 3 and 4 are all proportional to the logarithmic partner of this primary 
field (in a rank 2 Jordan cell). 

This should constitute a major step towards the exact solution of the Abelian sandpile model, 
since, if the conformal picture we propose is confirmed, it basically solves the aspects of the model 
concerned with the height variables, which are the most natural microscopic variables. (Indeed the 
simpler problem of identifying the scaling fields corresponding to the boundary height variables 
on open and closed boundaries was solved in [13, 14].) It however does not solve the part of the 
model which is related to other non-local variables, like avalanche observables, whose exact solution 
remains a challenge. 

Before getting there, we believe that a number of issues deserve further investigations, which 
lie in the context of the present article, and from which the sandpile model but also logarithmic 
conformal field theory in general could benefit. 

One of them concerns the unexpected nature of the logarithmic field describing the heights 2, 3 
and 4. We have argued that it has no realization in the free symplectic fermion field theory, widely 
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considered as the canonical local realization of a c = — 2 conformal theory. One certainly would 
like to know to what extent the sandpile realization and the symplectic fermion realization differ, 
and what the specific properties of the theory realized by the sandpile model are. 
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